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RESUMEN 

En este trabajo se estudian los posibles efectos locales que pueden presentarse en el 
emplazamiento donde se encuentra localizado el acelerografo de la ciudad de Olot 
(Girona). Como primer paso, se realiza una caracterizacion del subsuelo en el 
emplazamiento del acelerografo, basandose en estudios geotecnicos, pozos de agua y 
mapas geologicos existentes del lugar. A continuacion, se evalua la funcion de 
amplificacion del emplazamiento, analizandose la columna de suelo, mediante un metodo 
de modelizacion ID lineal-cquivalente, utilizando como senal de entrada un pulso de 
Ricker. La funcion obtenida se compara con el resultado obtenido en un estudio 
experimental realizado con registros de ruido ambiental. Finalmente se analiza la senal del 
sismo de Sant Pau de Fenolhet ocurrido en 1996 y registrada en el acelerografo situado en 
este emplazamiento. 



SUMMARY 

This work studies the possible local effects in the site of an accelerometer in the 
city of Olot (Gerona). A characterization of the subsoil of the site was performed based on 
geotechnical studies, lithology of water wells, and available geological maps. The transfer 
function is evaluated for the site using a soil column and modelizing it with an ID 
equivalent-linear method and a Ricker pulse as an input signal. The resulting transfer 
function is compared with results from a study of seismic noise. Finally, the signal from 
the 1996 Sant Pau de Fenolhet earthquake recorded in this site is analyzed using the soil 
column. 
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1 LOCALIZACION, GEOLOGIA Y GEOTECNIA DE LA CIUDAD OLOT 

La ciudad de Olot se encuentra localizada a 42° 10’ 54.8” Norte y 2° 29’ 36.3” Este, 
a una altitud de 443 m sobre el nivel del mar, en la comarca de la Garrotxa, provincia de 
Girona (Figura 1). La poblacion de Olot es de 27.482 hab., de acuerdo con el censo de 
1996 y tiene una superficie de 29.71 Km 2 (www.idescat.es). 

Olot es una zona densamente volcanica, con su consecuente variedad en las 
caracteristicas de la geologia del subsuelo. Las diversas erupciones volcanicas han 
propiciado el desarrollo de diferentes estratos de relleno, depositados en los lagos de la 
zona de esclusa volcanica (depositos lacustres y palustres). En muchos de los casos el 
tamano de estos depositos asociados al vulcanismo supera el espesor de los materiales 
volcanicos (Mallarach, 1 998). Las erupciones volcanicas van a crear condiciones idoneas 
para la acumulacion de sedimenlos al obstruir los valles y crear pequenas cuencas 
encerradas por conos o dentro de los crateres cerrados. Estos sedimentos lacustres y 
palustres se han intercalado con los materiales volcanicos y con otros sedimentos de origen 
fluvial o torrencial en depositos acumulados durante miles de anos. Por su extension, 
diversidad y riqueza de restos organicos estos sedimentos tienen un gran interes cientifico, 
ademas de su repercusion en el comportamiento dinamico del suelo. 




Figura 1 Localizacion de la Ciudad de Olot, las coordenadas utilizadas son las del sistema UTM. 



1.1 Informacion de la geologia v geotecnia de la ciudad de Olot 

Debido a las caracteristicas del subsuelo de Olot, fue necesario revisar la 
informacion disponible para poder caracterizar el subsuelo y en particular la columna del 
emplazamiento del acelerografo. De esta manera, la realizacion del estudio del subsuelo de 
Olot se hizo a partir de diferentes fuentes de informacion: mapas geologicos (Mallarach, 
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1982; ITGE, 1994), mapas topograficos (ICC, 1993), otros estudios (Fleta et al 1998), 
extraction de la base de datos del inventario de pozos de Catalunya (ACA, 2000) y 
estudios geotecnicos de la ciudad (Capdevila, 2000), Figuras 2 y 3. 




de estudio 



Figura 2. Localizacion del area de estudio. 




Kigura 3. Mapa geologico de la ciudad de Olot. 



Para esta primera parte se selecciono una area de 1.785 x 1.100 m (Figura 2), 
tomando como punto central la localizacion dc la estacion del acelerografo. Dentro de esta 
area se encuentran localizados 7 pozos de agua, y los 10 estudios geotecnicos disponibles 
(Figura 4). La information de los estudios geotecnicos proporcionan diversa information 
de la cual se pueden obtener las caracteristicas dinamicas de las diferentes capas del 
subsuelo de Olot, es ademas, mas completa que la proporcionada por la information de los 
pozos, sin embargo la profundidad maxima alcanzada es de 12.6 m, mientras que la 
profundidad maxima alcanzada en los pozos de agua es de 1 10 m. 




Figura 4. Area dc estudio de la ciudad de Olot para la caracterizacion del emplazamiento del acelerografo. 
Los niimeros en color rosa corresponden a los pozos de agua, mientras que los de color azul a los 
estudios geotecnicos. Se indican los trazos de los perfiles mostrados en la figuras 5 y 6. 
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De esta manera se realizo una interpretacion del subsuelo de Olot a partir de dos 
cortes geologicos indicados en la Figura 4 y representados en las Figuras 5 y 6, que se 
intersectan en la vertical del emplazamiento del acelerografo. En los cortes se pueden 
observar la presencia de varias coladas basalticas con potencias de varios metros, varias 
fallas subverticales y el reconocimiento de un aparato volcanico fosilizado no descrito 
anteriormente. La existencia de este aparato fosilizado podria explicarse por la intersection 
de dos fallas: una (SW-NE) que ya ha sido cartografiada y una segunda (NW-SE), no 
cartografiada cuyo trazado se apoya en la alineacion del volcan Montsacopa con los 
volcanes del Puig de la Costa y del Fontpobret. 

En la zona Ocste del rio Fluvia (ver figura 5) las primeras coladas, es decir las mas 
profundas, fueron producidas por el volcan la Garrinada (Mallarach, 1998). Por encima se 
tiene unas coladas posteriores, de las cuales no se sabe el origen. Segun Mallarach (1998), 
estas pcrtenecen a la Garrinada, pero por criterios geomorfologicos de estos edificios se 
deberia descartar esta ultima idea. Posterionnente ocurrio la formacion del edificio del 
volcan Montolivet y una colada producida por este que constituyen las primeras capas de 
P20 (ver figura 5). La formacion del volcan Montsacopa y la formacion de una colada con 
sentido hacia el Norte no afecta a los perfiles propuestos pero si los materiales piroclasticos 
que recubren el cono volcanico (parte superior sondeos S5 y S7). Cabe destacar en la parte 
mas occidental del perfil la presencia de basamento Eocenico. 




Figura 5. Perfil A, direction NW-SE 



En el perfil N-S mostrado en la figura 6 se pueden ver las mismas coladas 
mencionadas anteriormente y materiales piroclasticos del cono del volcan Montsacopa 
hacia el Norte.En resumcn, a pesar de la complejidad del subsuelo de Olot debido a la 
actividad volcanica, la columna del emplazamiento del acelerografo pucde defmirse con 
bastante fiabilidad hasta aproximadamente una profundidad de 50m, tal como se describe 
en el apartado siguiente. 
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520 N-S 



480 



440 



420 



Fjgura 6. Perfil B, direccion N-S. 




1.2 Columna obtenida a partir de la informacion geologica y geotecnica de Olot 

A partir de datos de los pozos de agua y de los estudios geotecnicos disponibles se 
ha definido una columna de suelo. La Tabla 1, muestra un resumen de las caracteristicas de 
la columna propuesta como son: el material, el espesor, la profundidad, el amortiguamiento 
(D %), el peso espedfico (ys) y la velocidad de las ondas de corte (Vs), asignados a cada 
uno de los estratos. Estas caracteristicas se obtuvieron de diferentes trabajos, entre ellos 
Cid, 1998, Zerva et al, 2000, Brown et al 2000 y Toshinawa et al , 2000, necesarias para el 
analisis con el programa lineal - equivalente SHAKE91 (Idris and Sun, 1992). 

La columna defmida para el emplazamiento del acelerografo consta de 6 capas con 
una profundidad maxima de 100 m, considerando a partir de ahi el basamento rocoso 
(Figura 8). La primera capa tiene un espesor de 2.5 m, compuesta con escorias y gredas 
como material principal. A esta capa se le asigno una velocidad de onda de corte de 150 
m/s. En los siguientes 35 metros se encontro la presencia de dos coladas basalticas. La 
primera colada es un basalto compacto con potencia de 21.5 metros, continuando con un 
basalto alterado o fisurado con una potencia de 3.0 m y por debajo de estas dos capas hasta 
llegar a los 38 metros se vuelve a tener basalto masivo. A los basaltos compactos se les 
asigno una velocidad de onda de corte de 900 m/s, mientras que al basalto alterado o 
fisurado se le asigno una velocidad de 550 m/s. La siguiente capa entre los 38 y 48 metros 
corresponde a un basalto fisurado, asignandole en este caso una velocidad de corte de 700 
m/s. La ultima capa corresponde a conglomerados, margas y gredas (capa 6), con una 
potencia de 52 metros (de acuerdo con ITGE, 94) y que llega hasta los 100 m de 
profundidad, a esta capa se le asigno una velocidad de 700 m/s. Todas estas capas se 
encuentran sobre un basamento rocoso, con velocidad de corte de 1 000 m/s. 
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Capa 


Material 


Espesor (m) 


D(%) 


ys (Ton/nv 1 ) 


Vs (m/s) 


1 


5 


2.50 


0.05 


1.43 


150.00 


2 


3 


21.50 


0.05 


2.40 


900.00 


3 


1 


3.00 


0.05 


2.30 


550.00 


4 


3 


11.00 


0.05 


2.40 


900.00 


5 


1 


10.00 


0.05 


2.30 


700.00 


6 


7 


52.00 


0.05 


2.30 


700.00 


7 


BASE 




0.05 


2.60 


1000.0 



Tabla 1. Litologia de la columna propuesta en el emplazamiento del acelerografo de la ciudad de Olot. 



(m) 

2.5(1 



3.00 



10.00 



BASE 




Escorias y gredas 




IVform-vion d* ctall. O) 

Basalto cnmpacto 




Basnlto fisurado 




DtfcnnarUnfe cUall* (%) 

Conglomerados y rnargas 



Figura 7. Columna de suelo utilizada correspondiente al emplazamiento del acelerografo junto con las 
curvas de variacion de paramentros dinamicos con la deformacion. 
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Para poder utilizar el programa SHAKE91, ademas de usar los datos anteriores es 
necesario definir las curvas dinamicas de los materiales que forman la columna de suelo. 
Estas curvas definen la variacion del modulo de cizalla dinamico (G) y del coeficiente de 
amortiguamiento (D) respecto a la deformation de cizalla (Figura 7 y 8). 

Las curvas utilizadas se escogieron de la tesis de Cid, 1998 y de las propuestas en el 
mismo programa SHAKE91. En la Figura 7, se muestra la columna, con los espesores de 
cada una de las capas, asi como sus correspondientes caracteristicas dinamicas de cada 
material. En la Figura 8, se muestran las curvas dinamicas de los materiales considerados, 
colocadas en una misma grafica para observar la diferencia entre ellas. Es importante 
senalar que existen curvas dinamicas en la literatura de una gran variedad de tipos de suelo, 
por lo que es posible utilizar 1 que mejor se aproxime al material considerado, cuando no 
se pueden realizar prucbas experimentales para la obtencion de ellas. 




Figura 8. Propiedades dinamicas usadas para los materiales de la columna propuesta indicadas en la Tabla 1. 



2 CALCULO DINAMICO DE LA RESPUESTA DEL SUELO DEL 

EMPLAZAMIENTO DEL ACELEROGRAFO DE OLOT 

La serial que se usara como entrada en el primer analisis es un pulso de Ricker con 
una amplitud maxima de 0.1 g, frecuencia predominante de 3 Hz y duration de 16 seg. Se 
ha considerado como valor maximo de la senal el de 0.1 g (Figura 9), de acuerdo al valor 
asignado a la zona de Olot en el mapa de peligrosidad sismica de Cataluna (Figura 10). 

Una vez que esta definida la columna de suelo, las caracteristicas dinamicas de los 
materiales y la serial de entrada, el analisis se realizara con el programa SHAKE91. El 
programa utiliza un metodo que considera la respuesta asociada a la propagation vertical 
de ondas de cizalla a traves del sistema viscoelastico lineal. La columna de suelo esta 
definida por N capas horizontales homogeneas e isotropicas de espesor A/, densidad 
modulo dinamico de cizalla G/ y coeficiente de amortiguamiento /?/. 
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Tletnpo (scr) 



Figira 9. Pulso dc Ricker con amplificacion maxima de 0.1 g. 




Figura 10. Mapa de peligrosidad sismica de Cataluna (Secanell et al , 1999). 

En este primer analisis se considera que la serial o pulso de Ricker, entra por la base 
de la col urn na litologica considerada como basamento rocoso (Outcrop) y se pretende 
calcular la amplificacion en la superficie. Las caracteristicas de los materiales utilizadas en 
el analisis se mostraron en el apartado anterior, sin embargo, de forma ilustrativa en la 
Figura 1 1 , se muestran las velocidades de corte para los estratos considerados y la forma en 
que van an con respecto a la profundidad. 



Vc Incidad ilc las ondas de corte (m/s eg) 

0 200 400 600 800 1000 




Figura 11. Velocidad de ondas de corte usadas para el analisis dinamico. 
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Una de las opciones que tiene el programa SHAKE91 es calcular el espectro de 
respuesta en cualquiera de la capas de la columna. En este caso se calcularon los espectros 
para la base de la columna o basamcnto rocoso y para la superficie. La Figura 12a, muestra 
el espectro de desplazamiento de ambos niveles, mientras que en la Figura 12b, se muestra 
los espectros de velocidad y finalmente en la Figura 12c, los espectros de aceleraciones. En 
estas tres graficos se puede observar claramente como se incrementan los valores de las 
respuestas entre el basamento y la superficie, indicando que existe una amplificacion de la 
senal,* como se explicara mas adelante. Los valores maximos para el espectro de 
desplazamiento se encuentra entre 0.5 y 1.0 Hz, mientras que en el espectro de velocidades 
se encuentran entre 0.8 y 1 .0 Hz, y entre 1 . 1 y 1 .4 los del espectro de aceleracion. 




Figura 12a. Comparacion de los espectros de Figura 12b. Comparacion de los espectros de 
desplazamiento entre las capas 1 y 7. velocidad entre las capas 1 y 7. 




Figura 12c. Comparacion de los espectros de aceleracion entre las capas 1 y 7. 



En la Figura 13, se ha representado las funciones de transferencia de cada capa con 
respecto al basamento rocoso. Es evidente que la amplificacion de la capa 7 con respecto a 
la misma es de 1, es decir no existe amplificacion alguna, de esta manera se comprueba que 
la senal ha sido considerada adecuadamente. Por otra parte se puede observar que la 
maxima amplificacion, en las frecuencias inferiores a 5 Hz, corresponde a la capa 6 con un 
valor de 1.77, esto es logico ya que las caracteristicas dinamicas de este material son 
menores que las del basalto masivo o fisurado y el espesor es de mas de 50 metros, esto 
provoca que la senal se amplifique. Sin embargo, la senal se reduce en las capas superiores, 
trabajando esta capa como un disipador de energia. La amplificacion maxima ocurre en la 
capa superficial, para frecuencias mayores de 10 Hz. 
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Finalmente, cn esta grafica se puedc observar que la amplification maxima en las 
frecuencias inferiores a 5 Hz, para la capa superficial es de 1.57, por lo que existira una 
amplification debil para estas frecuencias y que existe una amplification importante, 
superior a 5 solo para frecuencias entomo a 14 Hz. Asi pues la funcion de transferencia 
obtenida no muestra amplificaciones significativas en el rango de frecuencias de mas 
interes en ingenieria, es decir entre 0.1 y 10 Hz. 




Figura 13. Espectros de amplification de los estratos con respecto a la base de la columna. 

En la Figura 14, se muestran la senal de salida de la capa 6 y la capa 1, 
comparandolas con la senal de Ricker aplicada en el basamento. Se ha elegido la capa 6, 
por ser la que muestra la mayor amplification en los primeros 5 Hz. Si bien es cierto que la 
aceleracion maxima se presenta en la capa 1, la diferencia entre esta capa y la capa 6, es 
pequena, por lo que se puede comprobar que la amplification de las primeras 5 capas es 
mucho menor que la de la capa 6. 




Figura 14. Comparacion de la aceleracion entre la capa 1 (supcrficie) y la 7 (base). 
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Figura 15. Deformacion de los estratos. La grafica solo muestra el rango de 3.5 a 7 s, de 
la historia en el tiempo de la deformacion. 

Otra forma de observar el comportamiento de la columna del suelo, es revisando la 
deformacion que sufren las capas en la historia del tiempo. En la Figura 15, se puede ver 
que la capa que sufre la mayor deformacion es la 6, esto permite disipar energia y por lo 
tanto reducir la serial de salida en la superficie. 

Finalmente, en la Figura 16 se muestran a manera de resumen los valores maximos 
de esfuerzo, deformacion y aceleracion obtenidos en la columna. Los esfuerzos y las 
deformaciones son obtenidos a la mitad de cada capa, mientras que la aceleracion se 
obtienen en la parte superior de cada estrato. 



Esfuerzo de corte max (% ) Deformacion por corte max (% ) 

0 5 10 15 20 0 0.01 0.02 





Aceleracion maxima (g) 

0.09 0.1 Oil 0.12 0.13 




Figura 16. Valores maximos obtenidos en los diferentes estratos de la columna. 



12 



U. Mena, S. Figueras, J. Fleta, T. Susagna y .1. A. Canas 



3 COMPARACION CON LOS RESULT ADOS OBTENIDOS DEL ANALISIS 

EXPERIMENTAL A PARTIR DE REGISTROS DE RUIDO AMBIENTAL 

A parti r de registros de ruido con los cocientes espectrales de las componentes 
horizontales sobrc la vertical podemos analizar la posible existencia de un periodo 
predominante en la funcion de transferencia del suelo en el emplazamiento, es decir a este 
periodo apareccria una amplificacion del movimiento del suelo en caso de sismo. Asi podra 
comparase con los resultados obtenidos con los calculos realizados en la columna de suelo 
propuesta para el emplazamiento del acelerografo mostrada en el capitulo anterior. 

Para ello se ha realizado una experimentacion que consiste en medidas de ruido 
sismico en el emplazamiento de aproximadamente 3 minutos cada una. El equipo utilizado 
ha sido un acelerografo (K2 de la marca Kinemetrics ) de 19 bits de rango dinamico que 
permite una resolution de 40 pm/s 2 , suficiente para registrar los bajos niveles del 
movimiento del suelo correspondiente al ruido sismico. En los registros realizados los 
niveles alcanzados han sido, por termino medio, del orden de los 300 pm/s 2 . 

En la Figura 17 se muestran los analisis espectrales realizados. En la valoracion no 
se considerara los ruidos espureos que aparecen a frecuencias superiores a los 50 IIz que 
son debidos a las condiciones de la medida experimental (como por ejemplo la falta de 
suficiente aislamiento electrico del equipo). 

Se obscrva que para frecuencias inferiores a los 10 Hz la senal se mantiene 
practicamente liniforme sobre un valor del cociente de 1. Para frecuencias superiores 
aparece un pico en la serial a una frecuencia alrededor de los 14 Hz lo que indicaria una 
posible amplificacion a dicha frecuencia. 

Estos resultados estan de acuerdo con la funcion de transferencia (Figura 13) 
obtenida segiin la modelizacion en el capitulo anterior en la que no aparecen valores 
sign ificati vos de amplificacion para frecuencias inferiores a 10 Hz y solo aparece un pico 
entomo a 14 Hz, efecto de la capa delgada mas superficial. 
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4 ANALISIS DEL ACELEROGRAMA DE OLOT DEL SISMO DE SANT PAU 
DE FENOLHET DE 1996 

El dia 18 de febrero de 1996 se registro un terremoto a las 02 h 45 min (hora 
oficial) seguido de un gran numero de replicas, con epicentro en el Pirineo oriental ffances, 
siendo percibido en toda Cataluna (Butlleti Sismologic, 1996). La intensidad epicentral 
asignada de acuerdo a la escala MSK fue de VI. El sismo se registro en la estacion 
acelcrografica de la ciudad de Olot, con los valores maximos de aceleracion, velocidad y 
desplazamiento mostrados en la Tabla 2. 



Componente 


Aceleracion 


Velocidad 


Desplazamiento 


E 


3.83 cm/s 2 


0.004 g 


0.17 cm/s 


0.01 cm 


Z 


1.91 cm/s 2 


0.002 g 


0.13 cm/s 


0.02 cm 


S 


4.55 cm/s 2 


0.005 g 


0.23 cm/s 


0.02 cm 



l abia 2. Valores maximos de aceleracion, velocidad y desplazamiento de los Registros del 



sismo de Sant Pau de Fenolhet. 



El terremoto tuvo una aceleracion maxima de 4.55 cm/s 2 (0.00473 g), 
correspondiente a la componente horizontal N - S. En la Figura 18, se muestra el registro 
de la senal acotada entre los segundos 10 y 30 del registro total de 35 seg. Considerando 
que el intervalo de tiempo de muestreo para el registro de la senal fue de 0.005 s, se tiene 
un total de 4000 puntos que se utilizaran para el analisis de la columna de suelo. 




Figura 18. Registro de aceleracion en la componente S - N, de la senal de Sant Pau de 
Fenolhet. La senal se encuentra truncada entre los 10 - 30 s. 

Este registro obtenido en la superficie, se tomara como senal de entrada en la capa 
superficial y se utilizara la columna de suelo deducida en el capitulo 3, del primer analisis 
para obtener una senal en el basamento rocoso, es decir, se pretende deconvolucionar la 
senal registrada en la superficie, para obtener la senal que en teoria entraria por el 
basamento. En cuanto a las caracteristicas de los materiales y curvas dinamicas se 
consideraran las mismas que se utilizaron en el primer analisis. 
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Revisando los resultados se puede observer que son similares a los obtenidos con la 
senal de Ricker, es decir, los valores de los espectros en el basamento rocoso son menores 
en general que los oblenidos en la superficie. En las Figuras 19a, 19b y 19c, se muestran 
los espectros de respuesta de desplazamiento, velocidad y aceleracion respectivamenle (en 
funcion del periodo T). 

Comparando los espectros de desplazamiento se puede ver que en terminos 
gcnerales los valores en la capa superficial son mayores que los del basamento rocoso. En 
cuanlo a los valores de velocidad en cl rango de 0.25 a 0.40 segundos, los valores son 
mayores en el basamento, aunque en el resto de periodos los valores son mayores los de la 
capa superficial. Finaimente, en el espectro de aceleraciones existen dos picos uno en el 
periodo de 0.1 Is y el otro en el periodo de 0.29s. El valor maximo se obtiene en la capa 
superficial en el primer pico, aunque en el segundo pico el valor del basamento rocoso es 
mayor. El valor maximo del espectro de respuesta en superficie a un periodo 0.11s, es de 

0.018 g. 




Periodo (seg) 



Figura 19a. Comparacion de los espectros de 
deslazamiento entre la capa superficial (capa 
1) y el basamento rocoso (capa 7). 




Periodo (seg) 



Figura 19b. Comparacion de los espectros de 
velocidad entre la capa superficial (capal) y 
el basamento rocoso (capa 7). 




Figura 19c. Comparacion de los espectros de aceleracion entre la capa superficial (capa 1) 
y el basamento rocoso (capa 7). 
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En relacion a la funcion de transferencia que se puede observar en la Figura 20, se 
puede observar que son similares a los obtenidos en el primer analisis. La amplificacion 
maxima se obtuvo en la capa 6, correspondiente a los conglomerados y margas. Esta 
amplificacion fue de 1.77, a una frecuencia de 1.75 Hz, aunque la amplificacion obtenida 
en la capa 1, en esta misma frecuencia file de 1.57, por lo que se presenta el mismo 
fenomeno de desamplificacion de la senal obtenida en el analisis del primer ejemplo. 

El valor maximo obtenido para frecuencias mayores a 10 Hz, se alcanzo en la capa 
superficial con una amplificacion de 6.9 a una frecuencia de 1 5 Hz. Este valor es mayor 
que el obtenido en el ejemplo anterior que fue de 5.32, a una frecuencia de 14.38 Hz. 




Frecuencia (Hz) 



Figura 20. Espectros de amplificacion de los estratos con respecto a la base de la columna. 




Figura 21a. Historia en el tiempo de la aceleracion obtenida para las capas 1. 
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Figura 21b. Historia en el tiempo de la aceleracion obtenida para las capas 6. 




Figura 21c. Historia en el tiempo de la aceleracion obtenida para las capas 7. 



En la Figura 21a, 21b y 21c, se muestran las senates obtenidas en las capas 
seleccionadas de la columna de suelo teniendo en cuenta la funcion de transferencia 
deducida en el capitulo 3. Se puede observar que la serial obtenida en la capa 1, es la 
misma que se uso como serial de entrada. En la capa 6, la serial se reduce como 
consecucncia de la disipacion de energia. Y flnalmente se obtiene la senal en el basamento 
rocoso correspondiente a la capa 7. En esta capa el valor maximo obtenido fue de 0.00359 
g, menor que el de la senal de entrada de 0.00473 g. Esta senal corresponderia al registro 
en caso de tener el acelerografo en un emplazamiento rocoso. Como puede observarse, la 
modification es notable tanto por la reduccion del valor maximo de la aceleracion en un 25 
% como en la reduccion de contenido espectral de la senal por debajo de 0.2s, es decir, por 
encima de 8 Hz, como queda reflejado en los espectros de aceleracion de la figura 20c y en 
la propia senal (figuras 21 a y c). Cabe tambien sefialar que el pico de 0.01 8g observado en 
el espectro de aceleracion a 0.1 1 s (9 Hz) indica que la frecuencia de maxima amplification 
obtenida por calculo a 14 Hz podna ser mas baja, lo cual indicana un espesor superior de la 
capa de escorias bajo el emplazamiento. 
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Ticmpo (seg) 



Figura 22a. Historia en el tiempo de la deformacion en las capas 3. 




Figura 22b. Historia en el tiempo de la deformacion en las capas 6. 
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Figura 22c. Historia en el tiempo de la deformacion en las capas 7. 

En las Figuras 22a, 22b y 22c se muestran las deform aciones de las capas 3, 6 y 7. 
Se puede observar que las mayores deformaciones ocurren en la capa 6 (Figura 22b), por lo 
que ocurrc una mayor liberation de energia, reduciendo con esto la amplification de la 
serial . 
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Finalmente, en la Figura 23, se muestra los valores maximos de aceleracion 
obtenidos en cada una de las capas. Es posible ver en esta Figura, como el valor obtenido 
en el basamento rocoso es menor que el obtenido en la capa superficial. 



Aceleraci6n maxima (g) 
0.001 0.002 0.003 0.004 0.005 




Figura 23. Valores maximos obtenidos para las diferentes capas de la columna. 



5 CONCLUSIONES 

A parti r de la infomiacion geologica, de estudios geotecnicos y de las columnas 
litologicas de los pozos de agua, se ha caracterizado una columna de suelo en el 
emplazamiento del acelerografo. La profundidad maxima obtenida para esta columna fue 
de 100m hasta el basamento rocoso, dividida en 6 capas, con diferentes materiales entre 
ellos: escorias y gredas, basalto compacto, basalto fisurado, conglomerados y margas. En 
esta caracterizacion se constata la presencia de dos coladas de basalto compacto 
intercaladas entre conglomerados y basalto fisurado. 

La funcion de transferencia de la columna considerada no presenta amplificaciones 
significativas para las frecuencias inferiores a 10 Hz, solo en frecuencias mayores de 10 
Hz, se obtuvieron las mayores amplificaciones, debido a la pequena capa de escorias y 
gredas de 2.5 m. 

Los espectros estudiados para dos senales de entrada (Senal de Ricker y Registro de 
Sant Pau de Fenolhet), se pudo observar en terminos generales que son menores los 
obtenidos en el basamento rocoso que los obtenidos en la superficie, lo cual corrobora que 
efectivamente existe una amplification de la senal. 

Fue posible determinar la serial en el basamento rocoso, utilizando la columna de 
suelo definida para el emplazamiento del acelerografo e introduciendo el registro de Sant 
Pau de Fenolhet en la capa superficial. La aceleracion maxima en el basamento rocoso fue 
un 25% menor que el registro en superficie . 
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El espectro del registro en superficie indica una amplificacion importante entomo a 
9 Hz, ligeramente inferior a la frecuencia de 14 Hz obtenida por los otros metodos. Elio 
podria indicar que el espesor de la capa formada por escorias y gredas podria ser algo 
mayor. 



Finalmente, en la comparacion de los resultados determinados mediante el analisis 
numerico con los obtenidos en el analisis experimental con registros de ruido ambiente, se 
pudo comprobar, en ambos casos, que para frecuencias inferiores a 10 Hz, no aparecen 
valores significativos de amplificacion, mientras que los valores maximos se localizan en 
ambos casos entre las frecuencias de 1 0 y 20 Hz, debido al efecto de la capa delgada mas 
superficial. 
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RES l MEN 



En este trabajo se discuten algunos modelos que existen para representar el 
comportamiento aleatorio de las cargas vivas en las construcciones. Dichas cargas se suponen 
constituidas por cargas sostenidas mas cargas extraordinarias. 

Los modelos matematicos se representan mediante procesos estocasticos de Poisson. 
Por otro lado la magnitud de las cargas se representan por funciones de probabilidades tipo 
Gamma que se obtienen de muestreos hechos directamente en las construcciones. 

Usando los modelos matematicos descritos en la primera parte del escrito, se llega a 
valores de cargas maximas para el diseiio de oficinas y de salones de clase (escuelas). Estos 
valores contienen informacion sobre levantamientos de cargas sostenidas instantaneas 
realizados en inmuebles de la ciudad de Mexico; sin embargo, se considera que pueden 
aplicarse a construcciones de Latinoamerica, dado que el tipo de carga y los intervalos de 
tiempo en que se realizan cambios de ocupacion son similares a los de Mexico. 



1 INTRODUCE I ON 



En estudios previos (Ruiz y Aguilar, 1994) se llego a la conclusion de que la confiabilidad 
implicita en elementos de concreto reforzado disenados con el Reglamento de Construcciones 
para el Distrito Federal (RCDF) y con el reglamento del Instituto Americano del Concreto (ACI) era 
totalmente diferente. La confiabilidad se estimo a traves del indice p. En dichos estudios se 
demuestra que la confiabilidad de los elementos disenados con el reglamento mexicano es mucho 
mayor que la de los disenados con el reglamento americano. 

En dichos estudios se menciona que las grandes diferencias entre los valores de p 
asociados a uno y otro reglamento, se deben a distintas causas como son la diferencia en el valor 
de la resistencia del concreto f'c y la carga maxima de uno y otro reglamento (la carga adoptada 
por el ACI se estipula en el ANSI A58-1994). 
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Lo anterior motivo a la autora y a un grupo de estudiantes a estudiar con detalle las cargas 
maximas esperadas en diversos tipos de construcciones y los valores que se deben especificar en 
el RCDF. El presente constituye el primer trabajo realizado en Mexico con datos obtenidos 
directamente en oficinas y en salones de clase. 

En lo que sigue primeramente se describe la metodologfa que se siguio para realizar el 
estudio estadistico y posteriormente se aplica dicha metodologfa a oficinas y a salones de 
clase de varias escuelas publicas. Los resultados se usan para estimar los parametros de la 
funcion de distribucion de las cargas totales para diferentes areas y periodos de recurrencia de 
50 anos (oficinas) y 100 anos (escuelas). Los resultados se comparan con los que recomienda 
tanto el reglamento mexicano (RCDF) como el ANSI A58-1994 (que es la que adopta el 
reglamento del ACI). Para cada caso se dan recomendaciones sobre valores de cargas vivas 
maximas de diseno. 



2 REPRESENTACION DE LAS CARGAS VIVAS COMO PROCESOS 
ESTOCASTICOS 



Una estructura durante su vida util esta sometida a variaciones estocasticas de la carga 
viva. La magnitud y duracion de las cargas no se puede conocer en forma exacta, sin embargo 
es posible establecer modelos probabilistas que consideren su naturaleza aleatoria. 

El comportamiento de la carga viva en el tiempo que actuan en un area de piso se 
puede representar por dos tipos de cargas. Estas son: las sostenidas y las extraordinarias. A su 
vez cada una de estas se puede referir a cargas instantaneas o a cargas maximas. Enseguida 
se da una breve, explicacion de estas definiciones. Posteriormente se explica de manera 
general el proceso numerico que se utilizo en este estudio para calcular la distribucion de 
probabilidades de la carga maxima total. 



2.1 Carga sostenida instantanea 



Para el caso de las oficinas y los salones de clase las cargas vivas sostenidas 
instantaneas incluyen las debidas a papelerfa, mobiliario (pupitres, mesas, escritorios, etc.) y 
personas que laboran en las oficinas o bien estudiantes que asisten a clases. Su intensidad se 
considera constante durante ciertos intervalos de tiempo, hasta que ocurre un nuevo cambio de 
grupo de ocupantes. Se supone que estos cambios de ocupacion se pueden modelar mediante 
un proceso estocastico de Poisson. La variacion de las cargas en el tiempo se ilustra en la 
figura 1 .a. Se verified en este estudio que la funcion de distribucion de probabilidades que 
gobierna este tipo de cargas es tipo Gamma (Peir y Cornell, 1973). 

Los modelos de cargas vivas consideran que la intensidad de las cargas que se 
presentan en las estructuras actuan en forma estocastica en un area de piso. En este estudio 
se siguio el modelo propuesto por Ellinwood y Culver (1977). Estos autores proponen que la 
media es independiente del area y que la varianza no lo es. 

La hipotesis general es que las cargas sostenidas permanecen constantes durante 
ciertos periodos de tiempo hasta que exista un nuevo intervalo de ocupacion. Se supone que 
estos cambios de ocupacion tambien siguen un proceso de Poisson, por lo que el tiempo entre 
cambios de cargas sostenidas se puede modelar mediante una funcion de distribucion de 
probabilidades de tipo exponencial. 
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2.2 Carga sostenida maxima 



La distribucion de probabilidad de los maximos valores de la carga sostenida se adopto 
en este estudio como Extrema Tipo I (ver McGuire y Cornell, 1974). 



2.3 Carga extraordinaria instantanea 



Tanto en las oficinas como en las aulas de clase y en la mayoria de las estructuras, 
existen eventos que provocan una sobrecarga, estos se representan por las cargas 
extraordinarias, cuya magnitud puede alcanzar valores muy altos, pero su duracion (horas, dias 
o semanas) es pequena comparada con la de las cargas sostenidas, por lo que se considera 
que ocurren en un instante corto de tiempo (figura l.b). Este tipo de cargas se pueden 
presentar cuando existe la necesidad de llevar a cabo una remodelacion y se requiere 
almacenar el mobiliario y otros objetos en un area relativamente pequena con el fin de dejar 
espacio libre. Por ejemplo, se presentan incrementos considerables de carga en las oficinas y 
en los salones de clase cuando hay juntas o fiestas, es decir, existen aglomeraciones de 
personas en un area de piso. 

El modelo para cargas extraordinarias en este trabajo se representa a travds de 
eventos de ocurrencia independiente de Poisson. Durante cada evento se supone que las 
personas se reunen en grupos. El numero de cada uno de estos tambien puede representarse 
por un proceso de Poisson, con media que depende del area El modelo adoptado por McGuire 
y Cornell (1974) y por Ellinwood y Culver (1977) es el que se siguio en el presente trabajo. El 
valor esperado y la varianza son como sigue: 



E(E) = 



m w niR 



VarlE) = 



Xk(m R a\ 









En donde k convierte la carga a un valor de carga equivalente uniformemente distribuida y 
A es el area total. Los terminos m w , cr w , m R , <J R1 representan la media y la desviacion 

estandar, el subindice W se refiere al peso de una persona y R es el numero de personas en un 
grupo. 



CARGA 



J 



TIEMPO (AN<3S) 



Figura 1 . a Cargas sostenidas 
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Figura l.b Cargas extraordinarias 
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Figura l.c Cargas sostenidas mas extraordinarias 



Figura 1 . Representation de las cargas vivas. 



2.4 Carga extraordinaria maxima 

Se supone que la magnitud de las cargas extraordinarias se puede representar 
mediante una funcion de distribucion de probabilidades tipo Gamma. De igual manera que la 
carga sostenida maxima, los valores maximos de las cargas extraordinarias siguen una funcion 
de distribucion Extrema tipo I. 

El modelo de la carga extraordinaria maxima considera que el tiempo de espera entre 
cambios de carga sigue una distribucion exponencial, como consecuencia de suponer que el 
fenomeno esta regido por un proceso de Poisson. 

La funcion aproximada acumulativa (FL Xm (l)) de los valores maximos (Peir y Cornell, 1973) 
se expresa como: 



FUJI) = F e ( 1 ) e VJFMl,FMI 

En donde F E (I) es la funcion de distribucion de probabilidades Gamma acumulada de la 
carga extraordinaria instantanea. El parametro J e es la tasa promedio de ocurrencia de eventos 
extraordinarios y T el periodo de tiempo de interes (de acuerdo con la vida util de la estructura). 

Para niveles acumulativos altos, la ecuacion anterior puede aproximarse por (McGuire y 
Cornell, 1974): 
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FL Xm (I) 



e -v e T[l-F E (I )] 



Para obtener la media y la varianza de las cargas extraordinarias maximas en este 
estudio se siguio el modelo matematico propuesto por Y.K.Wen (1977). 



2.5 Carga total maxima 



La carga maxima originada por la combinacion de la carga sostenida y ia extraordinaria 
durante la vida util de !a estructura puede deberse a uno de los siguientes tres casos (segun 
Chalk y Corotis (1980). 

El primer caso considerado en este modelo es la suma de la carga maxima sostenida 
durante 1a vida util de la estructura (W s ) mas la carga extraordinaria mas grande que ocurre 
durante la carga maxima sostenida (W es ). 

Para el segundo caso se suman la carga maxima extraordinaria (W e ) que se presenta 
en la vida util de la estructura y la carga sostenida instantanea (W se ) que actua al mismo tiempo 
que W e . 



En adicion a estos dos casos se puede presentar un tercero. Este considera que se 
puede presentar, aunque con una menor probabilidad, la suma de la carga maxima 
extraordinaria durante la vida util de la estructura (W e ) y la carga sostenida maxima (W s ). 

La variacion de la carga total maxima debe considerar los tres casos expuestos 
anteriormente. En el presente trabajo se sigue el metodo de simulacion de Monte Carlo. Este 
considera de forma implicita estas tres combinaciones de carga, asf como la probabilidad de 
que se presente cada una de ellas. 



3 METODO DE SIMULACION DE MONTE CARLO 



El modelo de simulacion de cargas vivas por el metodo de Monte Carlo se resume en 
los siguientes pasos: 

1) Determinar los intervalos de tiempo entre cambios de carga sostenida. Los intervalos de 
duracion de las cargas sostenidas se generan a partir de una funcion de distribucion de 
probabilidades exponencial. Se generan intervalos de tiempo tj dentro del intervalo de tiempo 
de vida util (periodo de retorno T = 50 anos para oficinas, y 100 anos para aulas de escuelas), y 
un indice de cambios de ocupacion de cargas sostenidas (en este caso se adopto D s = 1 para 
salones de clase, y 1/8 para oficinas) . 

2) Generar la magnitud de la carga sostenida instantanea. Para generar valores de esta carga 
es necesario conocer la funcion de distribucion de este tipo de cargas. Esta se obtiene a partir 
de muestreos de las cargas instantaneas en los inmuebles. En estos casos la funcion de 
distribucion de la magnitud de las cargas instantaneas es tipo Gamma . 

3) Encontrar la carga extraordinaria maxima. Mediante la funcion de distribucion Extrema Tipo I 
y las formulas propuestas por Y. K. Wen (1979) se genera el valor medio y la desviacion 
estandar de la carga viva extraordinaria maxima. 

4) Elegir la combinacion mas desfavorable de la carga total. Esta es igual a suma de los 
valores de la carga sostenida instantanea mas la carga extraordinaria maxima (valores 
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generados en los incisos 2 y 3) durante los L intervalos de tiempo (valores generados en el 
inciso 1 ) que se presentan durante la vida util (T) de la estructura. 

5) Guardar la combinacion de carga total maxima mas desfavorable 

6) Repetir los pasos anteriores N veces. Donde N es el numero de simulaciones que se desea 
realizar (en este estudio N = 1000). 

7) Obtener los valores medios y desviaciones estandar de las cargas totales maximas 

8) Realizar los pasos 1 al 7 para distintas areas de influencia. 

El modelo de simulacion aqui expuesto obtiene de manera directa diferentes 
combinaciones de cargas vivas durante cada paso de simulacion. A partir de estas se obtiene la 
combinacion de valores maximos que pueden ocurrir en cierto intervalo de recurrencia. 



4 ESTUDIO ESTADISTICO DE LAS CARGAS 



En la primera fase de este estudio se obtuvieron histogramas de los objetos que mas 
frecuentemente se encuentran en las oficinas y salones de clase como son: libros, hojas de 
papel, folderes, sillas, computadoras personales, mesa bancos, etc. A partir de las mediciones 
hechas, se construyo una tabla basica de valores medios de aproximadamente 60 . objetos 
tipicos de oficina y muchos menos para salones de clase. Dicha tabla se modified dependiendo 
de las caracteristicas especificas de cada conjunto de inmuebies estudiados. 

El procedimiento general para la estimacion de cargas consistio de los siguientes 
pasos: 1) arreglos preliminares para solicitar autorizacion para visitar los inmuebies, 2) 
inspeccion inicial del inmueble y verificacion de su estado de carga, 3) peso de algunos 
muebles atfpicos, 4) eleccion del tamano de las areas basicas, 5) estimacion visual de las 
cargas instantaneas en las areas basicas, y 6) dibujo de la planta del inmueble en estudio. 

Los modelos probabilistas descritos anteriormente se alimentaron con datos 
estadisticos obtenidos directamente de oficinas y salones de clase en la ciudad de Mexico. Se 
muestreo un total de 14890 m 2 correspondiente a oficinas, y un area de 13700 m 2 a salones de 
clase. 

• El proceso anterior se aplico a los siguientes cinco edificios de oficinas: 

- Instituto de Ingenierfa, UNAM (Edif.2. Edificio mas cargado de los cinco tratados) 

- Direccion General de Inst. Tecnologicos (DGIT - DGETA) 

- Instituto Mexicano del Petroleo (SIPE y SIPPI) 

- Ingenieros Civiles Asociados (ICA) 



• Por otro lado, las escuelas visitadas fueron las siguientes: 

- Escuela Primaria “Insurgentes Bravo” 

- Escuela Primaria “General Francisco Menendez 

- Escuela Secundaria Diurna Num. 168 

- Escuela Secundaria Num. 185 

- Escuela Nacional Preparatoria Num. 5 “Jose Vasconcelos” 

- Escuela Nacional Preparatoria Num. 6 “Antonio Caso” 

- Facultad de Ingenierfa, Division de Ciencias Basicas, UNAM. 

- Facultad de Contadurfa y Administracion, UNAM. 
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4.1 RESULTADOS DEL AREA TOTAL ESTUDIADA CORRESPONDIENTE A 
OFICINAS 



Caraa viva sostenida instantanea. 



Los datos de las cargas sostenidas instantaneas de los cinco edificios de oficinas 
analizados se unieron y a partir de ellos se obtuvo la media (p) y la varianza (a 2 ) de la carga viva 
sostenida instantanea del area total estudiada (14890 m 2 ). 

= E[L] — 75.1 ^ ; 

m 



4647 74 

a 2 = Var[L]~ 409.29 + - 

A 

Se supone que los cambios de ocupacion de oficinas se rige por un proceso de Poisson, 
con tasa media de ocurrencia igual a 1/8 de ano. 




Caraa extraordinaria. 

En este estudio se consideraron los siguientes datos: (<J W , <7 W ) = (150, 25); W, 

representa el peso de una persona concentrada en una celda en lb, y (<7 R , cr R ) = (4, 2); R, 
representa el numero de personas cargando una celda. 

La ocurrencia.de las cargas extraordinarias tambien se supone como un proceso de 
Poisson, con tasa media de ocurrencia igual a 1 ano. 



Cargas totales. 

Con el fin de obtener la distribucion de las cargas totales maximas se apiico el modelo de 
simulacion de Monte Carlo descrito previamente en este trabajo. 

Los resultados de esperanza E y varianza (Var) de las cargas maximas totales se 
presentan en la Tabla 1 . La varianza muestra una reduccion en funcion del area. Esto ocurre tanto 
para cargas sostenidas instantaneas (s), como para cargas extraordinarias (e) y cargas totales 
maximas (t). 



En la figura 2 se presentan los valores que propone el ANSI-1994 y los valores medios 
que se obtienen con el programa de simulacion de Monte Carlo descrito previamente en este 
articulo. Como se ve en la figura 2, los resultados que se obtuvieron en este estudio resultan muy 
similares a los propuestos por el ANSI-1994. Los actuales valores de diseno especificados en el 
ANSI -1994 son muy cercanos a la media de la maxima combinacion de carga viva para el diseno 
de estructuras, asociada a un periodo de retorno de 50 anos. 

En la figura 3 se comparan estos valores con los propuestos en el actual Reglamento de 
Construcciones para el Distrito Federal (RCDF). La figura muestra en el eje de las abscisas el area 
tributaria. Debido a esto aparecen curvas relacionadas con vigas y con columnas. En dicha figura 
se puede observar que los valores del Reglamento RCDF resultan parecidos a los valores medios 
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que se obtienen de la simulation de Monte Carlo y a los valores nominales propuestos por el 
ANSI-1994 para areas tributarias pequenas; sin embargo, a medida que el area tributaria aumenta 
las diferencias son mayores. Es decir, los valores medios de las intensidades de cargas vivas 
totales que se obtuvieron en este estudio son significativamente menores a lo que propone el 
RCDF-1993, a pesar de que en esta investigation se revisaron cargas vivas en edificios de 
oficinas cuya densidad de ocupacion es alta. En la misma figura 3 se presentan los valores que se 
especificaban en el RCDF de 1976. Como se puede ver, estos valores son parecidos a los 
obtenidos en este estudio, para areas mayores que 36 m 2 . 

En la Tabla 2 se observa que los actuales valores nominales (W m ) del RCDF (1993) 
resultan muy conservadores, y que tienen probabilidades de excedencia P* cada vez mas 
pequenas a medida que aumenta el area de influencia, A|. Tambien se observa que los valores 
nominales (L) que propone el ANSI (1994) tienen probabilidades de excedencia (P*) cercanos a 
50% para areas relativamente pequenas, y probabilidades de excedencia ligeramente mayores 
(aunque no se alejan demasiado de 50%) a medida que aumenta el area. 



En la Tabla 2 a y u representan los parametros de la funcion Extrema Tipo I. La 
obtencion de estos parametros se describe mas adelante. 



TABLA 1 . Datos estadisticos y resultados de la simulacion. 



DATOS DE ENTRADA 



RESULTADOS DE LA SIMULACION 



AREA 

M 2 



18.6 

37.2 

74.3 
11T5 
148.6 
185.8 
223^0 



260.1 



kg/m 

75.1 

75~T 
75.1 
75. F 
75.1 
75.1 

75.1" 



75.1 



kg/m 

25.7 

23.1 
21 1 

21.2 
21 . 6 " 
20 . 8 " 



20.8 



20.7 



kg/m 
39 X 

45.7" 
37 7T 
31.4“ 
211 
25.1 



23.0 



21.4 



kg/m 

39.8" 

30.4 

19T 

14.6 

11.8 

~ioT 

8.8 



7.9 



E[e] 

kg/m 2 

159.5 

"1309 

87.1 

667" 

54.7 

47.2" 

41/7~ 



37.8 



Var[e] 

|(kg/m 2 ) 2 

1801.9" 

1043.4 

491 .9 

327.0 
249 1 

1 93 6~ 

149.1 

125.6 



EW 2 

kg/m 2 

90.5 

91 7 T 

~ 95^6” 
97.6"" 
99.3"" 
100 . 1 " 



100.8 



101.2 



Var[s] 

(kg/m 2 ) 2 

806.7 

652.1 

525.2 
469.9 
418/1 
385.5 
363.1 " 

351.4 



EP] 2 

kg/m 

250.0 

222 . 6 " 

182.2" 

164.3 
154 T 

147.4 



142.5 



138.9 



VarFtJ 

(kg/mT 

1734.5 

IOIOF" 

553. 9~" 

443.8 " 

403.2 

374.1 

353.6"" 

343.8"" 
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CARGA VIVA (kg/m 2 ) 




Area de infl. rm 2 i 


0.0 


18.6 


37.2 


74.3 


1115 


148.6 


185.8 


223.0 


260.1 


MOD. SIM. 5 EDIF. 




250.0 


222.0 


182.2 


164.3 


154.1 


147.4 


142.5 




ANSI-1994 — 


244.1 


244.1 


244.1 


190.5 


166.7 


152.6 


142.9 ' 


135.8 





Figura 2. Reduccion de las cargas vivas en funcion del area. 
(Datos de los cinco edificios analizados) 



CARGA VIVA (kg/m 2 ) 




AREA TR IBUTARI A (m 2 l 


o.o i 


9.3 


18.6 


37.2 


74.3 


1 1 1 .5 


148.6 


185.8 


223.0 


260.1 


RDF-93 — ' 


248.9 ! 


248.9 j 


248.9 


248.9 


228.7 


219.8 


214.4 


210.8 


208.1 


206.0 


RDF-76 


- 


i 257.8 


188.9 


168.7 


159.8 


154.4 


150.8 


148.1 


1 46.0 


ANSICOLUMNAS 


244.1 


244.1 1 


1 90.5 


152.6 


1 25.7 


1 1 3.9 


106.8 


101 .9 


98.4 


95.6 


ANSI VIGAS 




244.1 i 


244.1 | 


244.1 


1 90.5 


152.6 


1 35.7 


1 25.7 


1 18.9 


113.9 


1 10.0 


MOD ELO COL. 




222.2 j 


222.2 j 


182.6 


1 54.7 


1 34.9 


126.2 


120.9 


1 1 7.4 


: 114.7 


1 12.7 


MODELO VIGAS 


* 


222.2 ! 


222.2 1 


222.2 


182.6 


154.7 


142.3 


134.9 


129.9 


1 126.2 


1 23.3 



Figura 3. Reduccion de las cargas vivas en funcion del area. 
(Datos de los cinco edificios analizados) 
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TABLA 2. Probabilidades de excedencia para diferentes areas cargadas 





n 




a 




W m (RCDF) 

kg/m 2 


P 


■MjMfelM 


P 






1734.5 






250.0 


0.429 






37.2 










250.0 








74.3 








ygi 




0.013 






111.5 




443.8 


mi 
























0.004 


















0.0036 


142.9 


0.530 




142.5 


353.6 


0.0682 




219.8 


0.0029 


135.8 


0.588 




138.9 


343.8 


0.0691 


130.6 




0.0026 


130.2 


0.641 



De los parrafos anteriores se deduce que existe similitud entre los valores nominales del 
ANSI y los valores medios obtenidos en este estudio, y que estos discrepan de los especificados 
en el actual RCDF. 



5 ANALISIS DE UN EDIFICIO DE OFICINAS CON CARGAS MUY GRANDES 



En esta seccion se hace un analisis sobre las cargas del Edificio 2 del Instituto de 
Ingenieria de la UNAM. Se eiigio el estudio particular de este edificio por las siguientes dos 
razones: 1 ) la media de las cargas sostenidas instantaneas es mayor en este edificio que en todos 
los demas analizados, y 2) la variabilidad de las cargas con el area es menor que la que presentan 
los demas edificios que se estudiaron. 



CARGA VIVA (kg/m 2 ) 

300 -r- 




100 i * 



50 -} 



0 



AREA TRIBUITARIA (m 2 l 


0.0 


9.3 


18.6 


37.2 


74.3 


1 1 1.5 


148.6 


185.8 


223.0 ! 


260.1 


RDF - 1993 


- 


248.9 


248.9 


248.9 


248.9 


228.7 


219.8 


214.4 


210.8 


208.1 ! 


206.0 


MOD.SIM.5.EDIF.COLUMNAS — 


222.2 


222.2 


182.6 


154.7 


134.9 


126.2 


120.9 


117.4 


114.7 ! 


1 12.7 


MOD.SIM.5.EDIF.VIGAS 




222.2 


222.2 


222.2 


182.6 


154.7 


142.3 


134.9 


129.9 


126.2 | 


123.3 


MOD.SIM.ii. COLUMN AS 


4 #- 


246.9 


246.9 


209.6 


183.3 


164.7 


156.4 


151.5 


148.1 


145.7 ' 


143.7 


MOD.SIM.ii.VIGAS 


• 


246.9 


246.9 


246.9 


209.6 


183.3 


171.6 


164.7 


159.9 


156.4 


153.7 



Figura 4. Reduccion de las cargas vivas en funcion del area. 
(Datos de cinco edificios analizados, Instituto de Ingenieria y RCDF-1993) 
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Los resultados de este analisis se presentan en la figura 4. En esta se muestran varias 
curvas. Estas corresponden a lo siguiente: 1) los valores nominales especificados en el actual 
RCDF (Ifnea continua, sin sfmbolos); 2) los valores medios obtenidos con simulacion de Monte 
Carlo correspondientes a los datos de todos los edificios tratados en este estudio (MOD. SIM. 5 
EDIF); y 3) los valores medios obtenidos con simulacion pero considerando solo las plantas de 
Mecanica Aplicada y Estructuras del Instituto de Ingenieria (MOD. SIM. ii). La figura muestra que 
los valores medios de este ultimo caso se acercan a lo que propone el RCDF, sin embargo el 
RCDF presenta valores mayores. Es decir, que aun considerando areas de oficinas muy cargadas 
se obtienen valores menores a lo que establece el actual Reglamento. Lo anterior conduce a 
afirmar que las cargas especificadas en el RCDF son extremadamente conservadoras a medida 
que el area aumenta. 

Enseguida se encuentran leyes de reduccion de cargas totales con base en los resultados 
de este estudio. 



6 LEY DE REDUCCION DE LA CARGA VIVA CON EL AREA 



En lo que sigue se ajustan las curvas que gobiernan la reduccion de cargas con el area 
considerando los dos siguientes grupos de datos: 1) la estadfstica de los edificios en estudio (en la 
figura 5 se describe como: CINCO EDIFICIOS); y 2) la estadfstica solo del edificio mas cargado 
(llamado INSTITUTO DE INGENIERIA en la figura 5). Los datos del primer caso se refieren a las 
primeras dos columnas de la Tabla 2, excepto que en dicha tabla las areas son de influencia A ( , y 
en la figura 5 dichas areas son tributarias A T . El ajuste de las rectas se realizo utilizando la tecnica 
de mfnimos cuadrados. 



CARGA VIVA (kg/m 2 ) 




(AREA INFLUENCIA) A -y 2 



Figura 5. Formula de reduccion de cargas vivas en oficinas. 
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La recta ajustada con los datos correspondientes a todos los edificios en conjunto (se 
indica en la figura 5 como: CINCO EDIFICIOS) se muestra mediante una linea continua con 
pequenos cuadros. Su ecuacion es como sigue: 



L - Wm 



3.705'' 


/ 


0.392 + -j=- 

V yAr) 


V 



kg_ 

V m 2 



Donde Wm, es igual a 222.2 kg/m 2 ; A h es el area de influencia (en m 2 ). A[ es cuatro veces 
el area tributaria para una columna y dos veces el area tributaria para una viga. 

La recta que une las cruces corresponde a los datos del Instituto de Ingenierfa. Por otro 
lado, la recta continua, sin simbolos, se refiere a la ley de reduccion propuesta en este trabajo. Su 
ecuacion, valida para areas mayores que 36 m 2 , es la siguiente: 



donde A T , es el area tributaria 



W = 107 + 



852 
V At 



En la figura 6 se presentan las leyes de reduccion correspondientes a columnas y vigas de 
los cinco edificios analizados (Ifneas con cuadros muy pequenos, y con cruces); las asociadas al 
Instituto de Ingenierfa (con asteriscos, y cuadros un poco mas grandes); y la que se propone en 
este estudio (linea con rombos). En la misma figura se muestra la ley de reduccion especificada en 
el RCDF (linea sin simbolos). 



CARGA VIVA (kg/m 2 ) 




AREA TRIBUITARIA (m=l 


0.0 


9.3 


18.6 


37.2 


74.3 


111.5 


148.6 


185.8 


223.0 


260.1 


RDF - 1993 




248.9 


248.9 


248.9 


248.9 


228.7 


219.8 


214.4 


210.8 


208.1 


206.0 


MOD. SIM. 5. EDIF. COLUMNAS 




222.2 


222.2 


182.6 


154.7 


134.9 


126.2 


120.9 


117.4 


114.7 


112.7 


MOD.SIM.5.EDIF.VIGAS 


w 


222.2 


222.2 


222.2 


182.6 


154.7 


142.3 


134.9 


129.9 


126.2 


123.3 


MOD. SIM. ii. COLUMNAS 




246.9 


246.9 


209.6 


183.3 


164.7 


156.4 


151.5 


148.1 


145.7 


143.7 


MOD. SIM ii. VIGAS 


m 


246." 


246.9 


246.9 


209.6 


183.3 


171.6 


164.7 


159.9 


156.4 


153.7 


PROPUESTA 


♦ 


249 


249.8 


249.8 


249.8 


206.6 


188.5 


177.7 


170.3 


164.8 


160.6 



Figc a 6. Reduccion de las cargas vivas en funcion del area. 

(Datos de cinco edificios analizados, Instituto de Ingenierfa, RCDF-1993 y propuesta) 

En lo que sigue se da una breve descripcion de la ley de variacion que se propone en este 

trabajo. 
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Como se puede apreciar en la figura 6, para areas superiores a 36 m 2 , la carga es 
ligeramente superior a los valores correspondientes al analisis de superficies muy cargadas (como 
!a del Instituto de Ingenieria). Se propone que sea un poco superior a estos porque las cargas de 
dicho edificio podrian ser superadas. Esto se puede presentar ya sea porque se acumule 
demasiada papeleria de oficina en esta, o bien porque se cambie el uso del inmueble. Por ejemplo 
despues de los sismos de 1985 se descubrio que muchos edificios proyectados para oficinas se 
usaban como almacen de artfculos diferentes, como por ejemplo telas. 

En la zona cercana a los 36 m 2 la ley que se propone en este estudio, es conservadora. 
En esta zona se respeto lo propuesto en el RCDF de 1987: se conserva una carga uniforme hasta 
los 36 in 2 . Dentro de esta superficie se engioba la mayoria de los elementos secundarios y 
principales mas usuales en la practica. 

La carga uniforme para areas menores que 36 m 2 corresponde al valor medio obtenido 
por simulacion para 18.6m 2 (ver Tabla 2). En este caso es de 250 kg/m 2 . Es decir, coincide con lo 
especificado en el actual RCDF. 



7 CONCLUSIONES SOBRE CARGAS EN OFICINAS 



Las conclusiones de este estudio respecto a cargas vivas maximas en oficinas son las 
siguientes: 

a) Los resultados de esta investigacion muestran que los valores de carga nominal para oficinas 
que propone el actual Reglamento para las Construcciones del Distrito Federal (RCDF-1993) 
proporcionan una seguridad adecuada para areas pequenas, y que son muy conservadores para 
areas grandes. 

b) El RCDF recomienda valores nominales con probabilidades de^xcedencia muy bajas a medida 
que el area crece. Dichas probabilidades de excedencia, asociadas a un periodo de recurrencia de 
50 anos, varfan entre 42% (areas pequenas menores que 18.5 m 2 ), 0.70% (areas moderadas de 
aproximadamente 1 1 1 .5 m 2 ) y 0.26% (areas grandes de aproximadamente 260 m 2 ). Como se ve, 
estas probabilidades tienen variaciones muy bruscas con el area comparadas con las que se 
obtienen de los valores nominales especificados en el ANSI-1994, que son cercanos a 50% (ver 
Tabla 2). 

c) Comparando las cargas sostenidas medias de los edificios estudiados, se observan diferencias 
importantes entre un edificio y otro. En particular, las intensidades de las cargas sostenidas en una 
construccion (Edificio 2 del Instituto de Ingenieria) son significativamente mayores que las que 
corresponden a los otros casos. Procesando los datos que se obtuvieron para este caso se llego a 
expresiones para los valores de cargas medias totales en funcion del area, validas para superficies 
muy cargadas. Aun los valores medios obtenidos con estos datos son menores a los que se 
especifican en el actual RCDF. Mediante el estudio de este caso especial se corroboro que el 
RCDF es muy conservador para areas mayores que 36 m 2 aproximadamente. 

d) Segun lo que se observo despues de los sismos de 1985 y en este estudio, en ocasiones se 
cambia el uso del inmueble (parcial o totalmente), por lo que sus cargas actuantes no 
corresponden al tipo de ocupacion para el que fue disenado. Estas incertidumbres no se tomaron 
en cuenta directamente en los modelos probabilistas usados en este estudio; sin embargo, se 
trataron de considerar de manera muy simplificada mediante el analisis de un edificio muy cargado 
(Instituto de Ingenieria). 

e) Durante la realizacion de este estudio fue posible observarse de que en algunos edificios el 
espacio propio para oficinas se utiliza para bodegas, zona de archiveros, equipo de aire 
acondicionado, salas de computo y pequenas bibliotecas. Tomando esto en consideracion, se 
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propusieron valores de cargas medias totales en funcion del area superiores a los 
correspondientes a superficies muy cargadas. 

f) En este estudio se encontraron dos leyes de reduccion de la intensidad de las cargas vivas en 
funcion del area tributaria. 

1) La primera ley, basada en las estadisticas de cinco edificios analizados, asegura que 
las probabilidades de excedencia son aproximadamente del mismo orden (50%) tanto para areas 
pequenas como para grandes. 

2) La segunda ley toma en cuenta principalmente las estadisticas del edificio mas 
cargado. 

Esta ley es la que aquf se propone para diseno, en sustitucion de la que actualmente 
especifica el RCDF. Dicha ecuacion da lugar a valores mayores que los que se obtienen del 
analisis de un edificio muy cargado (Instituto de Ingenieria). Esto se considera asi con el fin de 
tomar en cuenta las incertidumbres debidas al cambio de uso del inmueble. 



8 RESULTADOS DE LA SIMULACION PARA EL CONJUNTO DE ESCUELAS 
ANALIZADAS 



Despues de realizar un analisis independiente para cada uno de los niveles educativos, 
se realizo un analisis tomando en cuenta el conjunto de escuelas. Para ello se unieron los 
datos de las cargas sostenidas instantaneas de las ocho escuelas anaiizadas se unieron. Con 
esos datos fue posible obtener la varianza en funcion del inverso del area. 

Los datos estadisticos para la carga viva sostenida instantanea de las ocho escuelas 
muestreadas son los siguientes: 



/i - 67.65 kg / rn l 
a 2 = 2307.8 / A (kg/m 2 ) 2 



En la tabla 3 se muestran los resultados de la simulacion para el conjunto de escuelas 
anaiizadas en forma simultanea. 

La forma de la distribucion de probabilidades (f.d.p.) de la carga total maxima se 
supone una funcion Extrema tipo I. 



En la figura 7 se muestra en forma comparative el valor medio de la carga maxima total 
en funcion del area para todas las escuelas en conjunto, para aquella que presento mayores 
valores medios (Escuelas de nivel superior) y los valores de diseno propuestos por la Sociedad 
Americana de Ingenieros Civiles (ASCE 7-95) y el Reglamento de Construcciones para el 
Distrito Federal (RCDF - 1993). 
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TABLA 3. Datos estadfsticos y resultados de la simulacion de todas las 
Escuelas analizadas 





DATOS DE ENTRADA 


RESULTADOS DE LA SIMULACION 


Area 

m 2 

(1) 


kg/m 2 

(2) 


°u 

kg/m 2 

(3) 


kg/m 2 

(4) 


kg/m 2 

(5) 


E[e] 

kg/m 2 

(6) 


Var[e] 

(kg/m 2 ) 

2 

(7) 


E[s] 

kg/m 2 

(8) 


Var[s] 

(kg/m 2 ) 

2 

(9) 


era 

kg/m 2 

(10) 


Var[t] 

(kg/m 2 ) 

2 

(11) 


37.16 


67.65 


7.897 


40.79 


27.16 


132.88 


665.55 


70.75 


73.42 


203.6 


614.68 


55.74 


67.65 


6.448 


36.64 


21.03 


121.38 


734.50 


70.19 


44.40 


191.5 


699.00 


74.32 


67.65 


5.582 


33.09 


17.30 


111.94 


779.91 


69.74 


32.43 


181.6 


737.85 


92.9 


67.65 


4.995 


30.30 


14.81 


104.08 


804.09 


69.19 


26.38 


173.2 


766.10 


111.48 


67.65 


4.559 


28.08 


13.01 


^96.997 


746.44 


68.92 


22.09 


165.9 


720.35 


130.06 


67.65 


4.222 


26.26 


11.65 


91 .728 


709.96 


68.73 


19.46 


160.4 


678.17 


148.64 


67.65 


3.948 


24.76 


10.58 


87.452 


707.16 


68.64 


16.68 


156.0 


675.10 


185.80 


67.65 


3.532 


22.39 


8.997 


79.761 


640.73 


68.58 


13.14 


148.3 


617.26 


222.96 


67.65 


3.224 


20.57 


7.877 


74.183 


624.57 


68.58 


11.29 


142.7 


596.48 


260.12 


67.65 


2.984 


19.15 


7.035 


68.737 


510.90 


68.45 


9.43 


137.1 


490.26 



9 PROBABILIDADES de excedencia para distintas areas 
CARGADAS 



Las probabilidades de excedencia se obtienen mediante la funcion distribucion Extrema 
TiP° ■. Ia °ual representa a la distribucion de la carga maxima total. La tabla 4 muestra en la 
segunda y terCera columnas los valores medios y la varianza de la carga total maxima para el 
conjunto de escuelas analizadas para distintas areas de influencia (indicada en la primera 
columna). Los valores de los parametros de la funcion de distribucion Extrema Tipo I se presentan 
en las columnas 4 y 5 de la misma tabla. Estos parametros son a y u, y se calculan con las 
siguientes ecuaciones: 



I 1.645 
y Var[jt] 



u = E\x\ 



0.577 

a 



Conocidos a y u, se puede determinar la probabilidad que tiene un cierto valor de 
carga de ser superado para un periodo de tiempo dado. 

A partir de los parametros y de fa forma de la distribucion antes mencionada se obtuvo 
la probabilidad de excedencia (P*) de los valores de carga maxima propuestos en el RCDF- 
1993 (W m ) y en las especificaciones de la Sociedad Americana de Ingenieros Civiles (ASCE 7- 
95). Dichas probabilidades se muestran en las columnas siete y nueve. 



36 



Sonia Ruiz y Alberto Soriano 




Figura 7. Valor medio de la carga maxima total del conjunto de escuelas.de las 
Facultades y valores de diseno de la ASCE 7-95 y del RCDF-1993 



Como se puede observar en la Tabla 4 la probabilidad de que se exceda el valor 
especifioado por el RCDF es muy pequena y decrece a 

probabilidad varia entre 0.0292 % (para areas menores que 400_ft ) y 0.0002 /« (para areas de 
2800 ft 2 ). Estas corresponden a un intervalo de interes de 100 anos. 

Por otro lado la probabilidad de excedencia de los valores especificados en el ASCE7- 
95 varian entre 57 % y 1.88% para las mismas condiciones antes mencionadas (ver columna 
l) E„ Se oaso ^ olL-va que sucede lo contrarlo, es decir las probabilidades de ex^deqcia 
son qrandes para areas pequenas aunque disminuyen para areas grandes. Para areas de 1000 
“ {qures fa^S toma iJL base la ASCE) la probabilidad de excedencia es de cas, 20 %, 
mientras que para el RCDF-1993 es de 0.0157 %. 

En la columna 7 de la tabla 5 se presentan probabilidades de excedencia para W m = 
250 kg/m 2 que es la carga maxima de diseno propuesta a partir de este estudio. 

En la figura 8 se presentan los valores especificados por el RCDF -1993 @50 ^m), 
los valores medios obtenidos para las escuelas de nivel superior los valores calculados p 
Harris Corotis y Bova (1981 ) y el valor propuesto er este trabajo (250 kg/m ). 
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TABLA4. Probabilidades de excedencia para diferentes areas cargadas 




W m 


P* 


L 


(kg/m 2 ) 




(kg/m 2 ) 


(RCDF 




(ASCE) 


(6) 


(7) 


(8) 


i 350 


0.00029 


195.6 


; 350 


0.00026 


195.6 



0.052 137.17 



0.053 131.78 




490.26 



0 I 0.00020 



0.00015 

0.00008 

0.00005 



0.00003 

ft , . 

0.00001 



0 0.00001 



350 0.00000 




195.6 


0.126 


195.6 


0.094 




195.6 




195.6 




195.6 


0.018 



TABLA 5. Probabilidades de excedencia para diferentes areas cargadas 
















































































38 



Sonia Ruiz y Alberto Soriano 




10 CONCLUSIONES Y RECOMENDACIONES PARA SALONES DE CLASE 



Con base en los resultados obtenidos del muestreo de cargas vivas en salones de 
clase y del modelo de simulacion Monte Carlo, se llega a las siguientes conclusiones respecto 
a las cargas vivas maximas en salones de clase: 

a) Los resultados de este estudio muestran que los valores de carga nominal que propone el 
Reglamento de Construcciones para el Distrito Federal (RCDF-1993) proporcionan una 
seguridad muy grande tanto para areas pequenas como para grandes. 

b) A partir de un analisis de los resultados obtenidos de la simulacion y que se muestran en la 
tabla 4 se pueden proponer valores de diseno, dependiendo de; 1) las probabilidades de 
excedencia que sea adecuado adoptar para que se tengan estructuras economicas y seguras, 
2) area de referencia que se seleccione, y 3) intervalo de tiempo de recurrencia que se 
considere adecuado. 

c) Las tablas 4 y 5 muestran las probabilidades de excedencia para la carga que proporciona 
el actual Reglamento mexicano (RCDF 1993), la carga recomendada por la Sociedad 
Americana de Ingenieros Civiles (ASCE 7-95) y la propuesta en este trabajo (250 kg/m 2 ). Sin 
embargo debe tenerse en cuenta que los factores de carga que se deben aplicar tienen que 
establecerse con las bases y niveles de confiabilidad adoptados en el Reglamento de 
Construcciones del Distrito Federal. Tambien debe considerarse que el modelo de cargas vivas 
empleado en el desarrollo del presente trabajo no considera efectos dinamicos provocados por 
cargas vibratorias sobre salones de clase. 

d) La carga viva maxima igual a 350 kg/m 2 (W m ) que recomienda el RCDF-1993, tiene una 
probabilidad de excedencia de 0.0292 % para un area de referencia de 37.16 m 2 (400 ft 2 ), 
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irnentras que una probabilidad de excedencia de 0.0157 % para un area de referenda de 92.9 
m (1000 ft ). Esta es el area que toma como referenda la Sociedad Americana de Ingenieros 
Civiles (ASCE 7-95) para proponer la carga viva de diseno (L = 195.6 kg/m 2 = 40 lb/ft 2 ). 

e) La carga propuesta igual a 250 kg/m 2 corresponde a una probabilidad de excedencia de 
2.2 % para un area de referencia de 74.32 m 2 , la cual corresponde al valor promedio de las 
aulas en el Distrito Federal, y a una probabilidad de excedencia de 1 .6 % para un area de 
referencia de 92.9 m (1 000 ft 2 ). 

f) En la figura 8 se muestran los valores especificados por el RCDF-1993 (350 kg/m 2 ), los 
valores medios obtenidos para las escuelas de nivel superior, el valor propuesto en la presente 
tesis (250 kg/m ) y los calculados por Harris, Corotis y Bova (1981). En esta figura se observa 
que los resultados obtenidos en el presente trabajo son muy parecidos a los de estos autores, 
aun cuando los datos corresponden a escuelas ubicadas en diferentes paises. La curva 
obtenida por Harris, Corotis y Bova (fig. 8) presenta un decremento mas pronunciado. 

g) Las probabilidades de excedencia para los valores especificados en el ASCE 7-95 y el 
RCDF-1993 considerando un periodo de recurrencia de 100 anos vari'an de 25.14 % para un 
area de 800 ft ( 74.32 m ) a 1.876 % para 2800 ft 2 (260.12 m 2 ), mientras que para el RCDF 
vari'an entre 0.02 % y 0.0002 %. 



11 CONCLUSIONES GENERALES 



Se han propuesto nuevos valores de cargas vivas maximas de diseno para el diseno de 
oficinas y de salones de clase. Ambos valores son menores a lo que actualmente se especifica en 
el reglamento mexicano RCDF-1993. Estos valores probablemente se modifiquen en futuras 
versiones del reglamento Mexicano. 

Es deseable analizar otras estructuras con diferentes tipos de ocupacion. Esto resulta 
conveniente para verificar la congruencia de la confiabilidad proporcionada en el RCDF por los 
valores de carga nominates para distintos usos de construcciones y areas tributarias. 
Actualmente se lleva a cabo un estudio similar al presente sobre estacionamientos. 
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RESUMEN 

En este trabajo se presenta un modelo que permite simular el comportamiento consti- 
tutive de materiales compuestos. El modelo se basa en la teori'a de mezclas basicas y se 
presenta una modificacion en la misma que permite tener en cuenta el fenomeno de des- 
lizamiento relativo entre las distintas fases que intervienen en el material. La verificacion 
del modelo constitutivo propuesto se realiza a traves de la simulacion numerica del com- 
portamiento de una probeta entallada de hormigon armado. En la simulacion numerica 
se presenta un estudio del fenomeno de transferencia de los esfuerzos entre la matriz (hor- 
migon) y el refuerzo (acero) y una comparacion entre un micromodelo, que presenta la 
discretizacion de las fases y un macromodelo, que utiliza el modelo constitutivo propuesto 
en este trabajo. 

ABSTRACT 

In this work a constitutive model for composite materials is presented. The model is 
based on the classical mixing theory and, due to the debonding phenomenon present in this 
type of materials, a modification is proposed. To check the constitutive model a numerical 
simulation of the behaviour of a reinforced concrete specimen with a notch in the central 
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-rva is presented. In the numerical simulation a study of the stress transfer phenomenon 
between matrix and reinforce phase is presented. A comparison between a macromodel and 
a micromodel of the specimen is given. 



1 INTRODUCCION 

En este artfculo se presenta un modelo matematico que permite simular el comportamiento 
constitutivo del hormigon armado. El modelo constitutivo esta basado en una teorfa general 
que permite tener en cuenta los diversos constituyentes que conforman un material. En 
particular para el caso de hormigon armado se consideran dos materiales: hormigon y 
acero. 

Esta teorfa general se denomina teorfa de mezclas de substancias basicas y permite 
tener en cuenta los diversos constituyentes de un material y la interaccion entre los mismos. 
La misma se basa, fundamentalmente, en la hipotesis de que la deformacion de todos los 
materiales que intervienen en el compuesto son iguales. Esta situacion es solo posible para 
el caso de materiales con comportamiento en paralelo, por. ejemplo barras de acero en una 
matriz de hormigon (hormigon armado), despreciando los fenomenos en los extremos de las 
barras de acero. Esta ecuacidjn de cierre, denominada tambien ecuacion de compatibilidad 
de deformaciones, no es valida para el caso en que se produzca un deslizamiento entre el acero 
y el hormigon. Este fenomeno se produce en el instante en que la tension cortante maxima 
admisible de la interfase entre el hormigon (matriz) y el acero (refuerzo) es superada. Por 
lo tanto es necesario introducir modificaciones en la teorfa general que permita tener en 
cuenta los deslizamientos relativos entre los constituyentes del hormigon armado. 

A continuacion se presenta en forma general los distintos modelos empleados para si- 
mular el comportamiento de materiales compuestos. Posteriormente, se introduce la teorfa 
de mezclas de substancias, que constituye el marco general que permite simular el com- 
portamiento de los materiales compuestos. En el apartado 4 se presenta la modificacion 
a la teorfa de mezclas que permite tener en cuenta el fenomeno de deslizamiento relativo 
entre las distintas fases del material. En el apartado 5 se presenta el modelo constitutivo 
elastoplastico utilizado en la matriz y para el zona de interfase se utiliza un modelo de 
daho de dano continuo (Oliver et ah, 1990) Por ultimo, en el apartado 6, se presenta una 
comparacion entre un micromodelo, en el cual se tienen en cuenta geometricamente cada 
uno de los constituyentes del material, y un macromodelo que utiliza la teorfa que aquf se 
presenta. 

2 MODELO CONSTITUTIVO PARA MATERIALES COM- 
PUESTOS 

Los materiales compuestos estan formados por diferentes tipos de sustancias inorganicas 
u organicas. En particular el hormigon armado se considera como un material compuesto 
bifase. 

Las caracterfsticas mecanicas de los materiales compuestos dependen de sus propiedades 
intrfnsecas y extrfnsecas como ser: caracterfsticas del proceso de fabricacion, tamaho de 
microporos y distribucion de los mismos, microfisuras, estados tensionales iniciales, etc. 
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Cada una de las sustancias componentes que integran el compuesto condicionan con 
su propia ley constitutiva el comportamiento del conjunto en funcion de la proporcion del 
volumen en que participan y de su distribucion morfologica en el compuesto. 

Existen diversas teorfas que permiten simular el comportamiento constitutivo de los 
materiales compuestos. La Teoria de Mezclas se considera adecuada para la simulacion del 
comportamiento de materiales compuestos en regimen lineal y con ciertas modificaciones 
permite simular, tambien, el comportamiento una vez superado el limite de proporcionalidad 
del material. 

La teoria de mezclas fue estudiada por Trusdell y Toupin (1960) y a su vez estos estudios 
dieron las bases de los trabajos de Ortiz y Popov (1982). Los resultados obtenidos por 
Trusdell constituyen tambien la base de trabajos como el de Green y Naghdi (1965) y el de 
Ortiz y Popov (1982) para materiales bifase. La teoria que aqui se presenta es mas general 
y representa el comportamiento constitutivo de un material compuesto de ”n” fases. 

La teoria de mezclas clasica generalizada para ”n” componentes permite considerar la 
combinacion simultanea del comportamiento de las diversas fases de un compuesto. Cada 
una de estas sustancias puede tener una ley constitutiva isotropa o anisotropa. La teoria de 
mezclas clasica es solo adecuada para simular ciertos materiales compuestos que presentan 
un comportamiento en paralelo, por ej. matrices con refuerzo de fibras largas (hormigon 
armado). A continuacion se presenta una breve introduccion a esta teoria. 

3 TEORIA DE MEZCLAS CLASICA 

La teoria de mezclas de sustancias basicas se basa en la mecanica del solido continuo local y 
se considera adecuada para explicar el comportamiento de un punto de un solido compuesto. 
Se basa en el prineipio de interaccion de sustancias componentes que constituyen el mate- 
rial, asumiendo las siguientes hipotesis: (i) en cada volumen infinitesimal de un compuesto 
participan un conjunto de sustancias componentes; (ii) cada componente contribuye en el 
comportamiento del compuesto en la misma proporcion que su participacion volumetrica; 
(iii) todos los componentes poseen la misma deformacion (ecuacion de cierre o compatibili- 
dad); (iv) el volumen ocupado por cada componente es mucho menor que el volumen total 
del compuesto. 

La segunda de las hipotesis implica una distribucion homogenea de todas las sustancias 
en una cierta region del compuesto. La interaccion entre las diferentes sustancias com- 
ponentes, cada una con su respectiva ley constitutiva, determina el comportamiento del 
material compuesto y depende basicamente del porcentaje en volumen ocupado por cada 
componente y de su distribucion en el compuesto. Esto permite combinar materiales con 
comportamientos diferenciados (elastico, elasto-plastico, etc.), donde cada uno de ellos pre- 
senta un comportamiento evolutivo gobernado por su propia ley. (Green and Naghdi, 1965) 
(Trusdell and Toupin, 1960) (Ortiz and Popov, 1982b) (Oiler et al., 1996). 

La tercera hipotesis establece que se cumple la siguiente condicion de compatibilidad 
bajo la hipotesis de pequenas deformaciones para cada una de las fases del material com- 
puesto: 



e ij"~ ( e ij) 1 — ( e ij ) 2 — — (£ij) n 



(1) 
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donde Eij y (s l j) n representan las deformaciones del conjunto y de la componente i-esima 
del material respectivamente. 

4 PROPUESTA PARA EL TRATAMIENTO DEL DESLI- 
ZAMIENTO ENTRE ACERO Y HORMIGON 

El objetivo del apartado es realizar modificaciones a la teori'a de mezclas presentada en el 
apartado anterior que permita tener en cuenta el deslizamiento entre las fases del hormigon 

armado (hormigon y acero). 

El analisis de los materiales compuestos en este trabajo se realiza utilizando la teorfa 
de mezclas de sustancias basicas. Esta teorfa se basa en la combinacion e inter accion de 
las sustancias basicas que conforman el compuesto. Se considera que en cada punto del 
material en analisis participan todas las sustancias componentes a la vez y cada una de 
ellas con su propia ley constitutiva y en la proporcion de volumen asignada. Todos lo 
componentes del material deben satisfacer la condicion de compatibilidad cinematica. Esta 
condicion establece que la deformacion total en un punto es igual para todos los componentes 
(Ecuacion de cierre). (ec. 1) 

Esta hipotesis es solo valida para materiales con comportamiento en paralelo, por ejem- 
plo matrices con fibras largas, despreciando los fenomenos en los extremos de las fibras. En 
el caso de materiales compuestos reforzados con fibras cortas la ecuacion de compatibilidad 
o cierre no se cumple. En este caso es necesario definir una ecuacion de cierre distinta (Oiler 
et al. ; 1995) o realizar una correccion en las propiedades de cada componente y mantener 
la ecuacion de cierre (Car et al., 1998). 

La ecuacion de compatibilidad tampoco es valida en el caso en que se produzca un 
deslizamiento relativo entre las fases del material (hormigon y acero en este caso). Este 
fenomeno se produce en el instante en que la tension cortante maxima admisible de la 
interfase entre acero y hormigon es superada. 

4.1 MODIFICACION DEL MODELO CONSTITUTIVO PARA CON- 
SIDERAR EL DESLIZAMIENTO ENTRE FASES 

Los materiales compuestos constituidos por una matriz .y un refuerzo poseen un comporta- 
miento no-lineal complejo, debido a que se produce un deslizamiento del refuerzo, motivado 
por la perdida de adherencia entre matriz y refuerzo. Este movimiento relativo entre re- 
fuerzo y matriz, produce tambien una perdida de rigidez en el conjunto, que se refleja en 
una disminucion de los parametros mecanicos del compuesto, sin que ocurra rotura en la 
fase de refuerzo (Hild et al., 1994) (Hild, 1994). 

La modificacion de la teorfa desarrollada en el apartado anterior, basada en la mecanica 
de medios continuos, consiste en introducir en la ecuacion constitutiva un comportamien- 
to basado fundamentalmente en el movimiento relativo de cuerpo rfgido que se verifica 
entre fibra y matriz. La incorporacion de este deslizamiento en la ecuacion constitutiva 
desarrollada en este trabajo se realiza teniendo en cuenta que este fenomeno presenta dos 
caracterfsticas: a) perdida de rigidez global debido a la disminucion de la colaboracion de 
1a, fibra en la matriz y b) deslizamiento relativo irrecuperable entre fibra y matriz. 
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Los materiales compuestos sometidos a estados tensionales en los cuales se ha produci- 
do el deslizamiento entre refuerzo y matriz no cumplen con la condicion de compatibilidad 
impuesta por la teoria de mezclas de sustancias basicas. Este fenomeno tiene como con- 
secuencia directa la imposibilidad de la matriz de transferir esfuerzos a la fase de refuerzo 
(acero). Esto es, el refuerzo no es capaz de aumentar su estado tensional por causas atri- 
buibles a la zona de interfase refuerzo- matriz. 

El modelo constitutivo se basa en la idea de que el proceso de transferencia de cargas 
de matriz a refuerzo varia en el momento en que la matriz sufre deformaciones plasticas. El 
movimiento relativo entre refuerzo y matriz puede representarse en mecanica de medios con- 
tinuos a traves de una deformacion plastica irrecuperable en el refuerzo. La determinacion 
del inicio de este fenomeno se realiza mediante una condicion de resistencia que compara 
la tension efectiva del punto con la resistencia de la fibra. Dada la forma en que participa 
el refuerzo dentro del compuesto y el mecanismo de transmision de tensiones entre refuerzo 
y matriz, la determinacion de su maxima resistencia o resistencia real y su capacidad de 
colaboracion depende de su propia resistencia nominal j ° resistencia del refuerzo 

en condiciones aisladas, de la resistencia nominal de la matriz (f N ) y de la resistencia 
nominal de la interface refuerzo- matriz (f N ) , o capacidad de transferencia de ten- 

V / ref —mat 

siones desde la matriz al refuerzo. Desde otro punto de vista, se puede decir que la fibra 
participa dentro del compuesto en funcion de su propia resistencia y de la capacidad de 
transferencia de esfuerzo de la interfaz fibra- matriz, por lo tanto su resistencia esta infiuen- 
ciada por el medio que la contiene. Se define entonces la resistencia del refuerzo contenido 
en una matriz como: 






/r 7 

t ref f V’’ /re/ J V* / mat V" / ref —mat 

donde Tf es el radio del refuerzo. La ec. (2) se deducen los siguientes casos limit es: 



( 2 ) 



• Si la matriz es mas resistente que el refuerzo y la adherencia refuerzo-matriz es perfec- 
ta, la capacidad de participacion del refuerzo queda limitada por su propia resistencia 
nominal (f*)^ = (/^ 

• Si se produce una falla en la matriz por microfisuras, etc., en tanto el refuerzo se 

mantiene en regimen lineal, la resistencia del refuerzo queda limitada por la resistencia 
de la matriz , pues no se podra transferir mas tension que la que le permite el medio 
que la contiene = (f N ') 

V / ref \ / mat 

• Si la falla se produce en la interface refuerzo-matriz, la resistencia del refuerzo queda 

limitada por la de la interface (f R ') = 2 (f N ^\ /^f 

\ / ref \ /ref— mat 

En la mayor la de los materiales compuestos se verifica que el agrietamiento por cortadura 
en la interfaz se produce antes que la rotura del refuerzo y se observa una separacion 
masiva entre refuerzo y matriz y por lo tanto la resistencia de la fibra queda limitada por 
la capacidad de la interfaz de transmitir esfuerzos. La aparicion de fenomenos plasticos en 
la matriz de un material compuesto sometido a un estado de cargas monotono creciente 




46 



E. Car S. Oiler E. Onate 



impide la transferencia de los esfuerzos desde la matriz hacia el refuerzo dando lugar a la 
aparicion de deformaciones irrecuperables por deslizamiento de la fase de refuerzo respecto 
de la matriz. A partir de este momento la transferencia de cargas de la matriz al refuerzo 
no es nula debido a la presencia de fenomenos de friccion entre ambas fases del material 
compuesto. Por lo tanto, el refuerzo aumenta su estado tensional segun un modulo elastico 
diferente del inicial. 

El modelo que se propone en este trabajo se basa en establecer el estado tensional 
en el que se encuentra el refuerzo en el momento de producirse la rotura en la matriz del 
compuesto. A partir de aqui las deformaciones irrecuperables debidas al deslizamiento entre 
ambas fases se tienen en cuenta con un modelo elastoplastico con endurecimiento para la 
fase de refuerzo debido a la presencia de las fuerzas de rozamiento entre fibra y matriz. 
Por lo tanto, se trata de un modelo no local material , debido a que el estado de uno de los 
materiales (refuerzo) depende de otro (matriz). 

A los fines de tener en cuenta la interaccion entre matriz y el refuerzo, es necesario 
determinar el estado tensional de la fase de refuerzo en el momento en que se verifica el 
deslizamiento relativo entre fibra y matriz debido a la aparicion de fenomenos plasticos en 
esta ultima. Este estado tensional es el maximo estado tensional que se transmite desde la 
matriz hacia el refuerzo bajo las condiciones de adherencia perfect a entre ambas fases. A 
partir de este momento los. esfuerzos se transmiten por fenomenos friccionales. En la Figura 
1 se muestra para el caso de un material compuesto formado por una matriz y fibra de 
refuerzos la interaccion entre ambas fases que permite determinar la resistencia del refuerzo 
en funcion de los criterios enunciados anteriormente. 

Con el objetivo de incluir el deslizamiento relativo entre fases en la formulacion del 
modelo constitutive para materiales compuestos se define la relacion 



r 



fr R 



( 3 ) 



donde representa la resistencia del refuerzo a n representa la tension en la direccion 
longitudinal del refuerzo en el momento en que se verifica el deslizamiento relativo entre 
refuerzo y matriz debido a la aparicion de fenomenos plasticos en esta ultima. 

A los fines de simular las deformaciones irrecuperables que se verifican como consecuencia 
del deslizamiento entre refuerzo y matriz es necesario redefinir el criterio de fluencia de la 
fase correspondiente al refuerzo. El factor r permite establecer el momento a partir del cual 
la matriz no es capaz de transferir carga al refuerzo y simular los fenomenos de friccion entre 
refuerzo y matriz. Este factor se considera constante a partir del momento de plastificacion 
de la matriz. 

Por lo expuesto anteriormente, para incluir el deslizamiento en el modelo constitutive 
es necesario redefinir el limite de proporcionalidad de las fibras, por lo tanto resulta: 



4>(t \g-,a) =F(tr ;g)-rk(a)= 0 (4) 

donde r e§ta dado por la ec. (3). El producto rk(a) permite redefinir el umbral de tensiones 
del refuerzo a partir del instante en el que se producen fenomenos plasticos en la matriz del 
material compuesto. La modificacion en la condicion de fluencia hace necesario redefinir 
nuevamente el modelo constitutive elastoplastico de la fase de refuerzo. 
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INCREMENTO: n | 




l 



PROXIMO 
INCREMENTO: n+1 



Figura 1: Esquema de solucion de problema no lineal bifase con DFM 
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5 EXPRESION DEL MODELO CONSTITUTIVO ELAS- 
TOPLASTICO DEL REFUERZO 

5.1 CONDICION DE FLUENCIA 

Segun se ha visto en el apartado anterior para tener en cuenta el deslizamiento relativo 
entre fases es necesario realizar una modificacion en la condicion de fluencia de la fase de 
refuerzo que permita redefinir la resistencia del refuerzo a partir del instante en el que se 
producen fenomenos plasticos en la matriz del material compuesto. 

El factor de proporcionalidad r afecta al limite de proporcionalidad del material. Este 
parametro no se introduce directamente en la condicion de fluencia con el objetivo de alterar 
lo menos posible la estructura matematica de la plasticidad La condicion de fluencia ec.(4) 
puede reescribirse como: 



4>{ct \ g\a) = l/(<r ;g) - K(a) = 0 (5) 

La ec.(5) representa el nuevo criterio de fluencia que permite tener en cuenta el fenomeno 
de deslizamiento entre refuerzo y matriz a partir de establecer el momento en el que la 
matriz no es capaz de transferir carga al refuerzo y simular los desplazamientos relativos 
irrecuperables entre refuerzo y matriz a traves de la plasticidad. 



5.2 REGLA DE FLUJO 

La regia de flujo estable la ley de evolucion de las deformaciones plasticas y se define como: 

L v {eY = a|£ (6) 

donde g = g{r\g) es la funcion de potencial plastico y A es una funcion escalar no negativa 
conocida como parametro de consistencia plastica que cumple con las condiciones de Kuhn- 
Tucker (Lubliner, 1990) (Crisfield, 1991). Para el caso de plasticidad asociada la ley de 
evolucion de las deformaciones plasticas queda: 



L v (e) p 



; 15 / 

r dr 



La teoria de la plasticidad exige tambien el cumplimiento de 



(7) 



A > 0 tj) (t ;g- a) < 0 A <j> (r \g\ a) = 0 (8) 

Las condiciones dadas por la ec.(8) se conocen en la terminologia clasica como condi- 

ciones de carga y descarga. Ademas ? debe cumplirse la condicion de persistencia expresada 
por la vari&cion temporal de la funcion de fluencia: 
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Material 1 

Matriz de Hormigon 


Material 2 
Refuerzo de Acero 


Material 3 
Junta Mat-R.efu 


Tipo de Comportamiento 


Mohr Coulomb 
Elasto-Plastico-Isbtropo 


Elastico-Isotropo 


Dano de Kachanov 


Mod. Young [kp/cm 2 ] 
Coef. Poisson 
Friccion interna 
Resist. Compr. [kpfcTTl 2 ] 
Resist. Tracci \kp/cm 2 } 
Gf , Gc . [kp/cm] 


3.5 x 10 5 
0.2 
30° 
200 
20 

0.25, 26.0 


2.1 x 10 e 
0.0 

2000 

2000 


3.5 x 10 5 
0.0 
30° 

20 

20 

2.0, 2.0 


Ley de comportamiento 
posterior a la fluencia 


Recto con 
ablandamicnto 




Exponencial con 
abiandamiento 



Tabla 1: Propiedades mecanicas de los materiales del micromodelo. 



6 EJEMPLO 

6.1 ESTUDIO DEL FENOMENO DE DESLIZAMIENTO ENTRE FA- 
SES A TRAVES DE UN MICROMODELO Y COMPARACION 
CON EL MODELO PROPUESTO 

En este apartado se presenta un ejemplo de aplicacion de la formulacion que combina la 
teorfa de mezclas, el modelo anisotropo en regimen de grandes deformaciones y la teorfa 
que permite tener en cuenta el fenomeno de deslizamiento fibra-matriz descrito en este 
trabajo. El ejemplo consiste en comparar la simulation numerica de una probeta de material 
compuesto (hormigon armado) con una entalla central sometida a un estado de traction en 
la que se han discretizado las fases de refuerzo y matriz (micromodelo), con una probeta 
similar en la que solo existe un material compuesto constituido por una fase de refuerzo y la 
matriz (macromodelo) . Las simulaciones se han realizado empleando una malla de elementos 
finitos rectangulares de 4 nodos con un total de 343 elementos, 392 nodos y 766 grados de 
libertad para el caso del micromodelo y 291 elementos, 336 nodos y 644 grados de libertad 
para el caso del macromodelo. En la Figura 2 se observan las mallas y las condiciones 
de borde empleadas para cada caso. En el micromodelo se considera que esta formado 
por tres materiales, matriz, zona de interfase fibra-matriz y refuerzo. El macromodelo 
esta constituido por un material compuesto constituido por dos fases: fibra de refuerzo y 
matriz. La Tabla 1 muestra las propiedades mecanicas de los materiales empleados en el 
micromodelo. Las propiedades mecanicas de las fases que constituyen el material compuesto 
del macromodelo son identicas a las de la matriz y refuerzo correspondientes al micromodelo. 

El objetivo del ejemplo es mostrar el fenomeno de transferencia de carga de la matriz 
hacia la fase de refuerzo. El estudio de este proceso se realiza a traves de un micromode- 
lo y luego se compara la respuesta carga-desplazamiento del mismo con un macromodelo 
constituido por un material compuesto y en el cual no se individualizan ffsicamente sus 
componentes. En la Figura 4 y 5 se muestra la evolution de la tension cortante en la zona 
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Figura 2: Malla de elementos finitos: a) micromodelo. b) macromodelo 



de interfase fibra-matriz para distintos incrementos de carga. En el primer incremento de 
carga, en el cual no se verifican procesos irreversibles, se observa que la distribucion de 
cortante a lo largo del refuerzo es similar a la curva teorica que se muestra en la Figura 3. 

En la zona central se observa un cambio en el signo de los esfuerzos debido fundamen- 
talmente a la presencia de la entalla. Las Figuras 6 y 7 muestran la evolucion de la tension 
longitudinal en la fase de refuerzo para distintos incrementos de carga. En la misma se 
observa, para el primer incremento de cargas, que en las zonas de los extremos del refuerzo 
se verifican las m&ximas tensiones cortantes mientras que las tensiones longitudinales crecen 
desde un valor nulo en el extremo hacia un valor constante a lo largo del refuerzo. En la zona 
central se observa, tambien, una variacion de la tension longitudinal debido a la presencia 
de la entalla. En estas figuras se observa, tambien, que un aumento en los esfuerzos aplica- 
dos produce fenomenos irreversibles en la zona de interfase .matriz-refuerzo en los extremos 
del refuerzo. Este proceso produce una disminucion de la capacidad de transferencia de 
esfuerzos desde la matriz hacia las fibras. Este fenomeno tambien provoca una modification 
en el estado tensional del refuerzo y se observa que la curva de distribucion de tensiones 
a lo largo del refuerzo deja de ser constante. En la Figura 8 se muestran las zonas de la 
interfase en las que se ha super ado el limit e de proporcionalidad del material para distint as 
etapas de carga. En la misma se observa que el fenomeno de deslizamiento relativo entre 
la matriz y el refuerzo comienza en la zona de los extremos de la fibra y progresa hacia el 
centro de la pieza. 

La Figura 9 muestra los contornos de desplazamientos en el primer y ultimo incremento 
de carga convergido para el micro y macromodelo. En el ultimo incremento de cargas se 
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Figura 3: Distribucion de esfuerzos axiales en la fibra y cort antes en la interfase fibra-matriz. 




•^Icremen. 1 
— — Icremen 10 

Icremen. 20 

—— Icremen. 30 
——Icremen. 40 
— Icremen 50 
Icremen. 60 
——Icremen. 70 
- Icremen. 80 
Icremen. 90 

“•Icremen, 

100 



Figura 4: Tensiones cortantes en la interfase fibra-matriz. Incremented 1-100 
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Figura 7 : Tensiones longitudinales en el refuerzo. Incrementos 110-300 

observa que los desplazamientos se concentran en la zona central de la probeta y a lo largo 
del refuerzo central. En los extremos de la pieza se puede apreciar el ” deslizamiento” entre 
fibra y refuerzo. 

La Figura 10 muestra las respuestas en fuerzas tot ales para ambos casos. Se observa que 
los resultados del micromodelo presentan un acuerdo satisfactorio con los del macromodelo. 
Es necesario destacar que el micromodelo no tiene la capacidad de simular los movimientos 
relativos entre las distintas fases pero tiene la ventaja de permitir simular refuerzos cuya 
discretizacion aumentarfa considerablemente el costo computacional del analisis debido a 
sus reducidas dimensiones respecto de las dimensiones de la pieza. 

7 CONCLUSIONES 

En este trabajo se presenta un modelo constitutive general para materiales compuestos 
basado en la teoria de mezclas de substancias basicas. Este tipo de modelos se basan en 
la hipotesis de que las deformaciones en todos los constituyentes del compuesto es igual y 
presentan el inconveniente de que son solo aplicables a materiales con comportamiento en 
paralelo bajo condiciones de perfecta adhesion. Existen materiales compuestos, bajo estados 
de carga superiores al limite de fluencia de la zona de interfase, en los que esta hipotesis de 
adhesion perfecta no se cumple gener£ndose, por lo tanto, desplazamientos relativos entre 
las fases del material. En el apartado 4 se presenta una propuesta de mejora que permite 
tener en cuenta este fenomeno. 

Por ultimo se presenta un ejemplo de comparacion entre un micromodelo y un ma- 
cromodelo de una pieza entallada de hormigon armado. El objetivo fundamental de esta 
comparacion es mostrar que a t raves del modelo constitutive que aquf se presenta, no es 
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Incremento 2 




Incremento 30 




Incremento 40 




Incremento 50 




Figura 8: Deformaciones plasticas en la interface fibra matriz para distintos incrementos de 
carga 
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Figura 10: Curva Fuerza - Desplazamiento macro y micromodelo. 



necesario realizar la discretizacion geometrica de cada uno de los materiales que constituyen 
el compuesto y constituye una potente herramienta que permite ahorro de coste computa- 
cional. 

Es necesario destacar que el micromodelo no posee la capacidad para simular los movi- 
mientos relativos entre las distintas fases pero tiene la ventaja de permitir simular refuerzos 
cuya discretizacion aumentaria considerablemente el costo computacional del analisis debido 
a sus reducidas dimensiones respecto de las dimensiones de la pieza. 
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RESUMEN 



Los paneles prefabricados de hormigon armado se usan en diversos lugares del mundo 
como sistema estructural para resistir acciones sismicas, debido a las indiscutibles ventajas que 
estos elementos presen tan desde el pun to de vista economico y de la reduccion de los plazos de 
construccion. Sin embargo, existe en la literatura escasa informacion sobre criterios de diseno. 
Los codigos sismorresistentes, en general, no incluyen especificaciones particulares para este 
tipo de construcciones. 

Se presentan en este trabajo recomendaciones practicas y un procedimiento propuesto 
para el diseno de paneles prefabricados, los que se vinculan entre si mediante un nuevo tipo de 
conexiones ductiles de acero. Dicho procedimiento se basa en los principios de diseno por 
capacidad y permite asegurar una adecuada disipacion de energfa, mediante deformaciones 
plasticas que se concentran en zonas especialmente disenadas, evitando la ocurrencia de otros 
modos de falla asociados con una respuesta fragil. Los paneles se conectan entre si con platinas 
de acero rectangulares, las cuales presentan una perforacion en el centro que asegura un 
comportamiento muy ductil de la union. Se incluyen ecuaciones simples para calcular la rigidez 
y resistencia de dichas conexiones, permitiendo al disenador plantear un modelo analftico 
simple de la estructura. 



1 INTRODUCCION 

El uso de paneles prefabricados de hormigon armado presenta significativas ventajas 
desde el punto de vista economico y de la reduccion de los plazos de ejecucion de obra. Es por 
ello que este tipo de estructuras se utiliza en todo el mundo para la construccion de edificios de 
baja y mediana altura. Cuando estos edificios se emplazan en zonas sismicas, tales como las 
zonas centro y oeste de la Argentina, se debe dar especial atencion al diseno de las uniones entre 
los distintos elementos que componen la estructura (fundacion, paneles, losas de piso y techo) 
con el objeto de evitar fallas de tipo fragil. La experiencia recogida en terremotos pasados 
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demuestra claramente que cuando este aspecto no es tenido en cuenta en forma adecuada, las 
consecuencias pueden ser muy graves, ocasionando dano severo, o aun colapso, en las 
estructuras con elementos prefabricados. Por ejemplo durante el terremoto de Northridge, 
California, 17/1/1994 [1] se observaron numerosos edificios construidos con paneles 
prefabricados con colapso parcial de la cubierta y colapso de los muros exteriores, dano que en 
la mayorfa de los casos se debio a problemas en las conexiones. 

Se presentan en este trabajo resultados experimentales y recomendaciones para el 
diseno sismorresistente de conexiones soldadas entre paneles adyacentes a los largo de las 
juntas verticales. Este estudio forma parte de un programa de investigacion general [2, 3] 
destinado a analizar el comportamiento sismico de edificios construidos con paneles 
prefabricados de hormigon armado. 



2 UNIONES ENTRE PANELES PREFABRICADOS EN LAS JUNTAS 
VERTICALES 

2.1 Generalidades 

Se emplean en la actualidad numerosos detalles diferentes para materializar las uniones 
entre paneles prefabricados a lo largo de las juntas verticales. Desde el punto de vista 
constructive, estas uniones pueden clasificarse en tres grupos: uniones abulonadas, uniones 
monolfticas con columnas adyacentes y uniones soldadas. Las Figs. A, B y C muestran algunos 
ejemplos tipicos de los tres grupos. 




(a) ( c ) 




Figura A. Uniones abulonados entre paneles prefabricados. 



Diseno sismorresistente para paneles prefabricados de hormigon armado con conexiones ductiles 



61 




(C) 



Figura B. Ejemplos usuales de conexiones monolfticas. 



Perfil anguto ^ Barras de refuerzo 

soldadas 




Soldadura 



Platina metalica 




(b) 



Figura C. Detalles de uniones soldadas. 
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Desde el punto de vista estructural, las uniones pueden dividirse en dos grupos. El 
primer grupo incluye aquellas conexiones que son capaces de trasmitir las fuerzas resultantes de 
la interaccion entre los paneles adyacentes, como es el caso de los detalles ilustrados en las Figs. 
B y C para conexiones monoliticas y soldadas. Por el contrario, el segundo grupo incluye las 
conexiones disenadas de modo que no se trasmitan fuerzas entre los paneles adyacentes, siendo 
su objetivo fundamental asegurar la estabilidad de los paneles prefabricados bajo acciones 
perpendiculares a su piano. Como ejemplo puede mencionarse el caso de uniones abulonadas 
en las cuales las platinas presentan agujeros sobredimensionados de forma tal que permiten el 
movimiento relativo de los paneles. 

La seleccion de la union mas adecuada para cada caso en particular depende no solo de 
consideraciones estructurales, sino tambien de las tecnicas constructivas empleadas y de las 
tolerancias admisibles en obra. 

Desde el punto de vista constructive, debe mencionarse que la practica usual es dejar un 
separacion de aproximadamente 10 a 15 mm entre los paneles adyacentes. Luego que los 
paneles has sido montados y conectados entre si, dicha junta debe sellarse a los efectos de 
asegurar la estanqueidad de la junta. Para lograr este objetivo se emplean normalmente distintos 
materiales disponibles comercialmente, los cuales deben reunir las caracterfsticas de ser 
elasticos (para absorber las deformaciones relativas) e incombustibles. 



2.2 Conexiones soldadas 

Las conexiones soldadas usuales estan formadas por platinas o perfiles de acero 
convenientemente. anclados en los paneles de hormigon, dispuestos de forma tal que puede 
soldarse a ellos una platina que actua de elemento de union entre ambos paneles (ver Fig. C). 
Los elementos criticos a considerar en una conexion de este tipo son: (i) la platina de conexion, 
(ii) las platinas o perfiles embebidos en el 'hormigon y sus anclajes, y (iii) los cordones de 
soldadura. 





Figura D. Platinas de acero usadas para conexiones soldadas ductiles. 



En construcciones sismorresistentes, las uniones soldadas deben disenarse de modo que 
exhiban una respuesta ductil, evitando los mecanismo de rotura fragil tales como los resultantes 
del inadecuado anclaje de las platinas de acero embebidas en el hormigon o de la fractura de los 
cordones de soldadura. Es conveniente que la fluencia (deformaciones plasticas) se concentre 
en una zona predeterminada de la conexion, la cual es especialmente disenada para tal fin. En la 
Fig. D se representan dos tipos de piatinas de acero propuestas para lograr este objetivo en los 
cuales se debilita una zona de la platina para crear un “fusible estructuraF’. La Fig. D (a) 
muestra el detalle propuesto en el programa de investigacion desarrollado por Restrepo et al. [2, 
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3], mientras que la Fig. D (b) describe el tipo de platina ensayado por Schultz et al. [4]. Ambos 
detalles son similares conceptualmente, sin embargo la primer platina descripta, Fig. D (a), 
presenta ventajas desde el punto de vista construed vo dado que es mas facil y economica de 
fabricar. 

Es importante destacar que en caso de que estas platinas fueran severamente danadas 
por la ocurrencia de un terremoto pueden reemplazarse facilmente, disminuyendo los costos de 
la reparacion y reduciendo el tiempo durante el cual la construccion no puede ocuparse. 

De acuerdo a los principios del diseno por capacidad [5], los elementos que pueden 
exhibir modos de falla fragil deben disenarse en funcion de la resistencia de fluencia de los 
elementos ductiles, considerando un factor de sobrerresistencia que tiene en cuenta un probable 
aumento de la resistencia real con respecto al valor nominal. En el caso particular de 
conexiones soldadas, la aplicacion de estos principios implica que los anclajes de las platinas o 
perfdes y los cordones de soldadura (elementos de pueden ocasionar una falla fragil) deben 
disenarse para resistir las acciones resultantes de la fluencia de la platina de conexion (elemento 
ductil). El factor de sobrerresistencia, X, depende fundamentalmente de las caractensticas del 
acero empleado para las platinas. A menos que se disponga informacion mas precisa se 
recomienda adoptar un valor de A.=1.25 para considerar un posible aumento de la tension de 
fluencia real en relacion al valor nominal considerado en el diseno, o bien determinar la 
resistencia ultima de la platina por metodos mas precisos. 




Figura E. Detalle del anclaje de las platinas de acero embebidas en los paneles. 



El anclaje de las platinas o perfiles de acero en el hormigon puede disenarse de modo 
que las fuerza de corte vertical, V s , resultante de la interaction de los paneles sea resistida por 
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dos barras de refuerzo colocadas en forma diagonal y soldadas a la platina o perfil, como se 
indica en la Fig. E. Es conveniente, ademas, disponer de dos barras horizontales para resistir 
cualquier momento adicional resultante de una posible excentricidad accidental. En base a estas 
hipotesis, el area A b de las barras diagonales debe ser como rmnimo: 



A b =X 



V.. 



2 / send 



( 1 ) 



donde f y es la tension de fluencia del acero de refuerzo y a es el angulo de inclinacion de las 
barras' diagonales. Ademas debe verificarse que los cordones de soldadura ce cada barra sean 
capaces de resistir la fuerza XA h f y , y que la longitud de las barra sea la necesaria para que se no 
produzca el deslizamiento por anclaje insuficiente. 

Los cordones de soldadura deben disenarse adecuadamente a los efectos de que tengan 
una longitud y espesor tal que aseguren que las platinas de conexion fluyen antes de que se 
produzca una fractura de los cordones, considerando los efectos de sobrerresistencia. 



3 RESPUESTA EXPERIMENTAL CICLICA DE PLATINAS DE ACERO 
PERFORADAS 



3.1 Introduccion 

El objetivo principal de este ensayo fue obtener informacion experimental sobre el 
comportamiento de platinas de acero con una perforacion circular en la parte central sometidas a 
esfuerzos de corte con variacion ciclica. Dicha informacion es de gran utilidad para verificar la 
respuesta del tipo de conexion propuesto (ver'Fig. D (a)) y para evaluar la rigidez y resistencia 
de la platina, dado que estos parametros son requeridos en el diseno de la union. 

Se presenta a continuation una breve descripcion del ensayo realizado y de los 
resultados obtenidos. Informacion mas detallada puede obtenerse del reporte final preparado al 
completar el programa de investigation [2]. 



3.2 Descripcion del ensayo 

Las caracteristicas geometricas de las platinas perforadas de acero (ver Fig. F) fueron 
adoptadas en base a un estudio analitico preliminar. Diferentes configuraciones fueron 
analizadas con un programa de elementos finitos en los que se aplicaron fuerzas de corte 
monotonicamente crecientes hasta producir la fluencia generalizada del acero. Los resultados 
del modelo numerico mostraron que se puede obtener una adecuada distribution de 
deformaciones plasticas alrededor de la perforacion circular si se adopta una relation entre el 
ancho de la platina de acero y el diametro del agujero de 1.6. 

La unidad de ensayo se construyo usando dos placas de acero con rigidizadores las 
cuales podian conectarse por medio de articulaciones a una maquina universal de ensayo con 
control de deformaciones. Las platinas de acero a ensayar, dos en este caso, se soldaron sobre 
cada una de las caras de las placas de acero con rigidizadores. Esta disposition se adopto para 
lograr una configuration simetrica y evitar movimientos fuera del piano. Las cargas axiales 
aplicadas en los extremos de la unidad de ensayo se transferian a traves de las platinas 
perforadas induciendo esfuerzos de corte en las misma. De esta forma se emulo el mecanismo 
de transferencia de corte entre paneles prefabricados sometidos a acciones sismicas. Una vista 
general de la unidad lista para ser ensayada se muestra en la Fig. G. 





El acero de las platinas perforadas se ensayo a traccion para determinar sus propiedades 
mecanicas. La tension de fluencia medida fue de 290 MPa y la tension ultima fue de 460 MPa, 
la cual se produjo con una deformacion de 23%. 

El ensayo se realizo aplicando las fuerzas en forma cuasi-estatica siguiendo una historia 
de desplazamientos ciclicos con amplitud creciente. Cada ciclo de desplazamiento se repitio 
dos veces hasta alcanzar una ductilidad de desplazamiento |i=12. A partir de este estado se 
aplico un desplazamiento creciente hasta |i=39 y se detuvo el ensayo, a pesar de que las platinas 
no mostraban degradacion de resistencia, dado que se excedio la capacidad de aplicar 
desplazamientos. 
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Debe mencionarse que las condiciones del ensayo no son totalmente iguales a las que 
experimentan las platinas de conexion en uniones reales. En primer lugar, las fuerzas fueron 
aplicadas en forma centrada, dado que se uso una configuracion simetrica con dos platinas 
perforadas, mientras que en construcciones con paneles prefabricados las platinas se colocan 
usualmente en un solo lado de los paneles. En segundo termino, las placas de acero sobre la que 
se soldaron las platinas a ensayar podian rotar libremente alrededor de dos articulaciones, a 
diferencia de lo que ocurre en conexiones reales donde las platinas de acero embebidas en los 
paneles prefabricados sufren rotaciones que dependen del momento flector actuante en los 
paneles y de sus caracteristicas. 



3.3 Resultados experimentales 

La comportamiento experimental de las platinas perforadas sometidas a corte demostro 
sin lugar a dudas que pueden proporcionar una respuesta altamente ductil cuando se las usa 
como elemento de conexion en uniones soldadas. La respuesta medida durante el ensayo, en 
terminos de fuerza de corte vs. desplazamiento vertical relativo, se presenta en la Fig. H. Puede 
observarse que los ciclos de histeresis son estables y presentan una area considerable, lo cual 
indica una buena capacidad para disipar energla. Ademas, las platinas fueron capaces de 
soportar desplazamientos relativos considerables (ductilidad de desplazamiento maxima p,=39) 
sin mostrar degradacion de rigidez o resistencia. 




Figura H. Respuesta experimental medida en el ensayo de carga ciclica. 



La fluencia del acero de las platinas pudo observarse a simple vista durante el ensayo. 
Las deformaciones plastic as se concentraron en la zona circundante al agujero central y se 
manifestaron en forma visible como line as de Liiders. La fluencia del acero comenzo en las 
zonas mas estrechas de las platinas debido a las deformaciones de flexion. A medida que se 
aumentaron los desplazamientos aplicados, la fluencia se propago hacia la zona central. En la 
Fig. I se presenta una fotografla de la unidad de ensayo al finalizar el mismo, pudiendo 
observarse que el orificio inicialmente circular se desformo significativamente, formando 
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aproximadamente una elipse. En el estado final, la relacion entre el eje mayor y menor de dicha 
elipse fue de 1.8, lo cual indica la magnitud de las deformaciones plasticas en esa zona. 




Figura I. Vista de las platinas de acero perforadas al finalizar el ensayo. 



El analisis de los cordones de soldadura realizado luego de la finalizacion del ensayo 
mostro que los mismos no presentaban dano aparente. 



3.4 Recomendaciones para el diseno de platinas de conexion ductiles 

Los resultados experimentales indican que las platinas con perforaciones ensayada en 
este programa de investigacion exhiben una respuesta ductil. La Fig. J presenta las dimensiones 
geometricas de la platina, longitud L, y ancho B, y de los cordones de soldadura en funcion del 
diametro de la perforacion. Se requiere realizar ensayos adicionales para verificar el 
comportamiento de platinas con diferente dimensiones relativas, pero hasta que esa information 
este disponible se recomienda utilizar el diseno propuesto en la Fig. J, cuyo comportamiento ya 
ha sido verificado experimentalmente. 

Consideraciones analiticas indican que el espesor de la platina, t, no afecta la 
distribution de tensiones en la platina, en aquellos casos en los que la hipotesis de tension plana 
es valida. Por esta razon, y tambien para limitar el espesor de los cordones de soldadura, se 
recomienda que la relacion ancho-espesor, B/t, se mantenga mayor que 5. 

El diseno de conexiones entre paneles de hormigon prefabricado requiere de la 
evaluation de la rigidez y resistencia de la union. Por ello se proponen ecuaciones desarrolladas 
en base a consideraciones analiticas y calibradas con los resultados experimentales, la cuales 
son de sencilla aplicacion y dependen de las caracteristicas mecanicas del acero y de la 
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geometria de la platina (los detalles de la deduccion de dichas ecuaciones pueden encontrarse en 
la referencia [2]), 




Figura J. Dimensiones relativas de la platina perforada de acero propuesta para las conexiones 

soldadas. 



La rigidez a corte de la platina perforada, K s , puede estimarse en base a un modelo 
simplificado que considera el efecto del agujero circular como una perforacion rectangular, 
mientras que la disminucion de rigidez debido a la flexibilidad de los cordones de soldadura se 
representa con una longitud equivalente de la platina mayor que la real. De esta forma puede 
demostrarse que: 

K s = 0.035 Et (2) 

donde E es el modulo de elasticidad del acero. 

Para la evaluation de la fuerza de corte, V sy , que produce la fluencia de la platina se 
considera que la fluencia se inicia en las zonas cercanas a la perforacion debido a las 
deformaciones de flexion, de modo que dichas zonas pueden considerarse como vigas en 
voladizo con seccion variable. En base a esta hipotesis se demuestra que: 

V sy =022tDf y (3) 

La resistencia a corte ultima depende del aumento de tensiones en el material debido a 
los efectos de endurecimiento de post-fluencia y a la propagation de las deformaciones plastica 
hacia la zona central de la platina. Asumiendo una distribution simplificada de las tensiones de 
corte en la seccion central de la platina, se encuentra que: 

V m =0.2StDf u (4) 

donde f u es la tension ultima o de rotura del acero. 

Es muy importante aclarar que si se emplea la Ec. 4 para calcular la resistencia de la 
platina, el factor de sobrerresistencia X debe tomarse como 1.0., dado que dicha ecuacion 
representa la resistencia ultima de la platina. 
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4 RESPUESTA EXPERIMENTAL DE DOS PANELES PREFABRICADOS 
CONECTADOS CON PLATINAS PERFORADAS 



4.1 Introduction 

Actualmente se utilizan una gran variedad de uniones diferentes para conectar paneles 
de hormigon armado en construcciones prefabricadas. Sin embargo existe escasa information 
experimental sobre el comportamientos de esas uniones bajo la accion de cargas sismicas. 
Debido al uso extendido de conexiones soldadas, se realizo un ensayo de una estructura 
completa formada por dos paneles prefabricados, los cuales se vincularon entre si con las 
platinas perforadas descriptas en la section previa. Ademas debe destacarse que este tipo de 
conexion exhibe un comportamiento relativamente complejo (comparado, por ejemplo con 
uniones abulonadas), dado que deben trasmitir la fuerza de corte resultante de la interaction de 
los paneles. 



4.2 Description del ensayo 

La unidad de ensayo se diseno para representar parte de la estructura tipica de un 
edificio de dos pisos, con una altura total de 7.0 m, con una losa prefabricada en el primer nivel 
y una cubierta metalica liviana. Para facilitar la construccion de la estructura y la aplicacion de 
cargas, la unidad de ensayo se realizo a escala reducida 1:2. La Fig. K ilustra las caracteristicas 
principales de la unidad de ensayo y sus dimensiones. Como puede observarse en dicha figura, 
la estructura esta formada por dos paneles prefabricados montados sobre un base, la cual se fijo 
a la losa del laboratorio mediante pemos de acero de alta resistencia para simular la condition 
de empotramiento. La losa del primer nivel y la cubierta metalica no fueron modeladas en el 
ensayo. Los paneles, de 75 mm de espesor, 2.0 m de longitud y 3.5 m de altura, se construyeron 
con hormigon, cuya resistencia a compresion al momento del ensayo fue de 28.8 MPa. En el 
centra de los paneles se coloco una malla de acero electrosoldado, formada por barras de 4 mm 
de diametro con una separation de 75 mm. El acero de dichas barras fue ensayado a traction 
arrojando una resistencia de fluencia de 415 MPa. Ademas se dispusieron armaduras especiales 
de refuerzo en las esquinas inferiores de los paneles, donde se materializa la conexion con la 
base. 



Los paneles se conectaron entre si con dos uniones soldadas con platinas perforadas 
(125 x 40 x 5 mm), cuyas dimensiones relativas se adoptaron segun el criterio indicado en la 
Fig. J. Las platinas de acero embebidas en el hormigon durante la construccion de los paneles 
se construyeron siguiendo el criterio ilustrado en la Fig. E. Las uniones se dispusieron en una 
cara de los paneles. 

La conexion entre los paneles y la base se realizo por medio de barras de refuerzo de 8 
mm de diametro, convenientemente ancladas en la base y que sobresalen en forma vertical, las 
cuales Se insertan en ductos metalicos dispuestos en los paneles de hormigon armado. Una vez 
que los paneles son montados sobre la base, los ductos se rellenan con “grout”. Este tipo de 
union tambien fue estudiada experimentalmente en > un ensayo preliminar (los resultados, 
criterios de diseno y conclusiones de dicho ensayo pueden encontrarse en las referencias [2, 3]). 

La Fig. L ilustra el portico de carga utilizado para aplicar las acciones laterales que 
simularon la accion sismica. Puede observase, ademas, que se emplearon dos vigas de acero, 
ubicadas a los costados de los paneles, para evitar movimientos fuera del piano de los mismos. 
Desplazamientos y deformaciones en puntos caracteristicos de la estructura fueron medidos con 
potenciometros lineales y bandas extensometricas, mientras que las acciones aplicadas por los 
actuadores hidraulicos fueron registradas con celdas de carga. Todos los instrumentos de 
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medicion se conectaron a un sistema de adquisicion de datos que permitio almacenar y anaiizar 
la informacion en forma digital en una computadora personal. 




Figura K. Caracteristicas principals y dimens i ones de la unidad de ensayo. 




Figura L. Portico de carga usado para el ensayo de los paneles acoplados. 



Las fuerzas laterales se aplicaron en forma cuasi-estatica siguiendo una historia ciclica 
de desplazamientos. Cada ciclo, con amplitud creciente, se repitio tres veces para imponer 
condiciones de carga severa que mostraran claramente la presencia de efectos tales como 
degradacion de rigidez y resistencia. El nivel de aplicacion de las fuerzas laterales se determino 
considerando la resultante de las fuerzas inerciales originadas sobre los paneles, losa del primer 
nivel y cubierta. 
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Para compensar los efectos de la reduction de escala en la construction del modelo (que 
altera la relation entre area transversal y volumen del panel) se agregaron 92 lingotes de plomo, 
los cuales se montaron sobre perfiles angulo colocados sobre cada cara de los paneles. 

4.3 Evaluation de la respuesta 

Resulta de gran interes determinar la respuesta esperada de la estructura a ensayar, para 
lo cual se desarrollo un modelo analitico basado en consideraciones simple de diseno (ver Fig. 
M). A1 aplicar una fuerza lateral V, a una altura h, cada panel prefabricado queda sometido a un 
momento flector en la base igual a: 

M=Vh-V s L (5) 

donde V s es la fuerza resistida por cada una de las platinas de conexion y L es la longitud de 
cada panel. 




Figura M. Acciones sobre los paneles prefabricados. 

La carga axial actuante en cada tabique se ve afectada por la interaction de los paneles 
entre sf, la cual origina que las platinas transfieran la fuerza de corte V s (desde el punto de vista 
conceptual, la estructura es similar a dos tabiques de hormigon armado acoplados por vigas). 
Por lo tanto la fuerza axial N en los paneles 1 y 2 es: 



N x = N g - 2V S 
N 2 = N g +2V s 



( 6 ) 




72 



Francisco Crisafulli, Jose Restrepo y Robert Park 



siendo N g la carga axial originada por las acciones gravitatorias en cada panel. 

La resistencia a flexion de las conexiones panel-fundacion puede determinarse mediante 
analisis sectional, teniendo en cuenta las caracterfsticas de los materiales, las position y tipo de 
barras de refuerzo y el esfuerzo axial en cada panel. En el caso considerado, cada panel tiene 
una resistencia diferente, M\ y M* 2 , debido a que las cargas axiales netas son diferentes (ver 
Ecuacion 6). Una vez que se han determinando los momentos resistentes se puede evaluar la 
resistencia lateral de la estructura, V*, asumiendo que es igual a la suma de las resistencias 
individuales de cada panel, de modo que: 

* M*\ + M*2 +2V L 

V — \i) 

h 

Para el caso de la estructura ensayada, V*=71.1kN. 

El procedimiento descripto previamente representa un solution simple para evaluar la 
resistencia de la estructura basado en consideraciones de equilibrio. El desplazamiento de 
fluencia y ultimo, tambien puede evaluarse mediante un procedimiento simplificado, 
considerando los desplazamientos originados por la deformacion elastica de los paneles mas los 
efectos de la rotation resultante de la deformacion en las barras de refuerzo que atraviesan la 
conexion entre el panel y su base de fundacion [2]. 

Ademas, se realizo un estudio mas detallado del la respuesta de los paneles utilizando 
un modelo de elementos finitos. Cada panel se represento con elementos de cuatro nodos, en 
tension plana, asumiendo comportamiento lineal y elastico del material. La junta vertical entre 
los paneles y las juntas horizontales entre cada panel y la base de fundacion se modelo con 
elementos tipo interfaces de cuatro nodos. Para representar el comportamiento axial y de corte 
(accion de dovela) de las armaduras de acero que cruzan la conexion panel-fundacion se 
emplearon resortes no lineales. Este mismo tipo de elementos se usaron para modelar las 
platinas de conexion entre paneles. 



4.4 Resultados 

En el ciclo de carga inicial, los paneles conectados con platinas ductiles se comportaron 
en rango lineal, sin mostrar fisuras o dano alguno. En el segundo ciclo, con una fuerza lateral 
maxima de 60 kN, se desarrollaron fisuras a lo largo de la junta horizontal entre los paneles y la 
base. Mediciones realizadas durante el ensayo muestran que las barras de acero en la conexion 
panel-fundacion fluyeron en el segundo ciclo de carga, al igual que las platinas perforadas de 
acero que conectaban los dos paneles. A medida que se continuo el ensayo aumentando los 
desplazamientos laterales se observo que los paneles sufrian deslizamiento por corte en la junta 
horizontal, sin embargo dicho deslizamiento fue controlado por las barras de refuerzo que 
cruzaban la conexion panel-fundacion. El ensayo se finalizo por la fractura de una de las barras 
de acero que conectaban uno de los paneles a la base de fundacion. La inspection visual 
realizada demostro que no habia fisuras en el panel prefabricado ni dano alguno en la zona de 
las uniones entre paneles a lo largo de la junta vertical. 

Los resultados del ensayo indican que la estructura exhibio un respuesta estable bajo la 
accion de cargas laterales con variation ticlica. En el rango inelastico, los ciclos de histeresis 
mostraron un estrechamiento significativo (“pinching”) como resultado del deslizamiento por 
corte de los paneles prefabricados. A pesar de este comportamiento, la estructura fue capaz de 
soportar las acciones laterales sin t degradation de resistencia, aun cuando se aplicaron 
desplazamientos laterales de importancia (la ductilidad maxima impuesta fue |i=22.7).. La Fig. 
N presenta la curva fuerza lateral-desplazamiento medida durante el ensayo. 
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Figura N. Respuesta fuerza lateral-desplazamiento medida durante el ensayo de dos paneles 

con uniones soldadas ductiles. 



En la Fig. N se presenta una comparacion entre los resultados experimentales con la 
prediction simplificada pudiendose observar, en general, una buena concordancia. La curva de 
resultados numericos calculada con el modelo refinado (la cual no se incluye en la Fig, N) fue 
capaz de representar en forma mucho mas precisa la respuesta experimental. Sin embargo, debe 
aclararse que una modelacion analitica con elementos finitos en campo no lineal no se utiliza 
usualmente en diseno debido a su complejidad. 

El desplazamiento lateral que experimentaron los paneles prefabricados puede 
desglosarse en tres componentes: (i) desplazamiento debido a deformaciones de flexion y corte, 
(ii) desplazamiento debido a deslizamiento por corte en la junta horizontal, y (iii) 
desplazamiento resultante de la rotacion rfgida del panel originada por la elongation de las 
armaduras de refuerzo vertical en la conexion panel -fundacion. El analisis de la information 
recogida durante el ensayo permitio estimar las distintas componentes, llegandose a la 
conclusion que en los ciclos iniciales la mayor parte del desplazamiento lateral fue resultado de 
la rotacion del panel, mientras que en los ciclos finales el deslizamiento por corte fue el 
mecanismo mas importante. 

Los instrumentos de medicion que registraron el desplazamiento relativo entre los dos 
paneles permitio evaluar la curva fuerza-desplazamiento para cada platina perforada de acero. 
Estos resultados fueron totalmente similares a los obtenidos en el ensayo de las platinas aisladas 
(ver Fig. H), lo cual confirma la gran capacidad de este tipo de conexiones para deformarse en 
forma ductil y estable. 



5 DISENO DE MUROS PREFABRICADOS CON CONEXIONES DUCTILES 

Los paneles de hormigon prefabricado con conexiones ductiles soldadas se comportan, 
desde un punto de vista estructural, en forma similar a tabiques de hormigon armado acoplados. 
El acoplamiento se logra a traves de las platinas de conexion, lo cual origina una diminution 
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del momento flector en los paneles y modifica la fuerza axial actuante en los mismos. La 
magnitud de estos cambios depende del grado de acoplamiento, esto es de la fuerza de corte V s 
que las platinas pueden transferir. A pesar de que ciertos componentes del sistema estructural 
pueden exhibir un comportamiento fragil, la aplicacion de los principios fundamentales del 
diseno por capacidad asegura que pueda obtenerse una respuesta global de la estructura 
adecuada para resistir acciones sismicas. 

A los efectos determinar las fuerzas sismicas actuantes en la estructura, y de acuerdo a 
los criterios presentes en las mayoria de los codigos sismorresistentes en vigencia, el disenador 
debe adoptar un valor de ductilidad global de la construccion. Teniendo en cuenta las resultados 
experimentales obtenidos en este trabajo, en los que se observo una respuesta aceptablemente 
ductil, pero con un mecanismo controlado por deslizamiento por corte, se recomienda disenar 
este tipo de construcciones como estructuras de ductilidad moderada, esto es con un factor de 
ductilidad global no mayor de 3.0. 

Se presenta a continuacion, en forma resumida, los pasos sucesivos del procedimiento 
de diseno propuesto para edificios con paneles de hormigon prefabricado con conexiones 
soldadas ductiles: 

1. Los paneles prefabricados deben disenarse para resistir los distintos estados de carga 
resultante de las fuerzas gravitatorias, de acciones laterales (en el piano del panel y en la 
direction perpendicular) y de esfuerzos termicos. Ademas, se debe prestar especial 
atencion, a los estados de carga resultantes del transporte y montaje de los paneles 
prefabricados. 

2. La resistencia a flexion de la conexion entre cada panel prefabricado y la base o fundacion 
debe analizarse considerando el efecto de la fuerza de corte V s transmitida por la union 
soldada de acuerdo a la configuration particular de cada caso (las Ecuaciones 5 y 6 
representan el caso de la estructura ensayada). Inicialmente, la fuerza V s debe ser estimada 
por el disenador con el criterio de evitar que la modification de la fuerza axial en los 
paneles conduzca a situaciones inconvenientes, tales como aumento o reduction excesiva de 
la resistencia flexional. Se recomienda evitar la situation de que alguno de los paneles 
prefabricados quede sometido esfuerzos de traction. 

3. El esfuerzo de corte resistido por cada panel, mayorado por los probables efectos de 
sobrerresistencia, tiene que ser transferido a la base o viga de fundacion. Para ello debe 
verificarse que la conexion panel-fundacion presenta la resistencia a corte adecuada. El 
procedimiento de diseno del panel y su conexion con la base (pasos 1, 2 y 3) se describe con 
mayor detalle en las referencias [2, 3], 

4. Las platinas perforadas de acero y los cordones de soldadura usadas en las conexiones 
deben seleccionarse siguiendo el criterio descripto en este trabajo (ver Fig. J), utilizando las 
Ecuaciones 2, 3 y 4 para calcular la rigidez y resistencia a corte de las mismas. En esta 
etapa deben seleccionarse las platinas de conexion para que tengan una resistencia similar a 
la adoptada en el paso 2. Caso contrario sera necesario repetir el proceso de diseno desde el 
paso 2. 

5. Por ultimo, se debe verificar que las platinas de acero embebidas en los paneles (sobre las 
que se sueldan las platinas perforadas) sean capaces de resistir las fuerzas de corte que se 
desarrollan en las platinas de conexion. Para ello se recomienda el uso del detalle de 
anclaje ilustrado en la Fig. E, siendo necesario determinar el area de las barras diagonales 
mediante la Ecuacion 1. Ademas deben verificarse los cordones de soldadura y la longitud 
de anclaje de las barras. 
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6 CONCLUSIONES 

Las principales conclusiones obtenidas en el estudio analitico-experimental que se 

presentan en este trabajo son: 

♦ Es de fundamental importancia disenar adecuadamente las uniones entre los distintos 
componentes de una construccion con paneles prefabricados de hormigon armado. Las 
conexiones representan los elementos mas debiles de la cadena de resistencia y deben 
disponer de sufic iente reserva de ductilidad para resistir acciones sismicas sin ocasionar el 
colapso de la estructura o danos irreparables. 

♦ Para el caso particular de conexiones soldadas entre paneles se propone en este trabajo el 
uso de platinas perforadas de acero. Estudios analiticos basados en modelos con elementos 
finitos y verificaciones experimentales indican que este tipo de elementos de conexion 
permiten obtener una respuesta sumamente ductil y estable, con buena disipacion de 
energfa. Para ello es necesario evitar la ocurrencia de modos de falla fragil, tales como los 
ocasionados por fractura de los cordones de soldadura o falla de los anclajes de las platinas 
embebidas en el hormigon. Este objetivo puede lograrse mediante la aplicacion de los 
principios basicos del diseno por capacidad. 

♦ A los efectos de verificar el comportamiento de la conexion propuesta, se ensayo la platina 
perforada de acero en forma aislada sometida a esfuerzos de corte con variacion ciclica. 
Ademas, se realizo el ensayo de una estructura completa formada por dos paneles 
prefabricados, montados sobre una base o viga de fundacion, conectados entre si con 
platinas perforadas de acero. 

♦ Se presenta un procedimiento simple para el diseno sismorresistente de paneles 
prefabricados con conexiones ductiles, basado en consideraciones analiticas y en los 
resultados obtenidos del programa experimental. La aplicacion de los principios basicos del 
diseno por capacidad asegura que, a pesar de que algunos componentes del sistema 
estructural puede exhibir un comportamiento fragil, la respuesta global de la estructura 
adecuadamente ductil. 

♦ Las platinas de conexion propuestas son faciles de fabricar y pueden reemplazarse en forma 
rapida y economica en caso de que fueran danadas con resultado de la ocurrencia de un 
sismo severo. 

♦ Mas investigacion experimental es necesaria para verificar el comportamiento de otra 
uniones comunmente usadas en la construccion de edificios con paneles prefabricados. 
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NOTACION 

B : ancho de la platina de conexion 
D : diametro del agujero circular en la platina de conexion 
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f y : tension de fluencia del acero 

f u : tension ultima del acero 

K s : rigidez a corte de la platina perforada 

h : altura de aplicacion de la fuerza lateral 

L : longitud de la platina de conexion 

L w : longitud de los cordones de soldadura horizontales 

M ' : momento Hector 

M* : memento resistente 

N : carga axial 

N g : carga axial originada por las acciones gravitatorias 

t : espesor de la platina de conexion 

t w : espesor de los cordones de soldadura 

V* : resistencia lateral 

V s : fuerza de corte en la platina de conexion 

V sy : resistencia de fluencia de la platina de conexion 

V su : resistencia ultima de la platina de conexion 

a : angulo de inclinacion de las barras de anclaje diagonales 

X : factor de sobrerresistencia 

p : ductilidad de desplazamiento 
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RESUMEN 

En este trabajo se presenta un modelo rapido y sencillo de evaluacion de la 
vulnerabilidad sismica de puentes a escala regional. El modelo fue estructurado de forma 
que puede ser aplicado especialmente a zonas donde no se cuenta con information de danos 
sfsmicos de puentes. El metodo se basa en la identification de las caracterfsticas mas 
relevantes e influyentes en el dano que sufrira un puente bajo la action del sismo. La 
valoracion de estas caractensticas se realiza mediante la determination de diecinueve (19) 
parametros, a los cuales se les asigno un grado de vulnerabilidad y un valor de importancia 
relativa en base a la opinion de expertos. Debido a las imprecisiones y subjetividad de las 
opiniones, estas fueron procesadas mediante tecnicas difusas. Una vez calificados cada uno 
de los parametros, se precede a relacionar las calificaciones de los parametros con su valor 
de importancia por medio de un promedio ponderado difuso. El resultado de esta relation 
por medio de tecnicas basadas en conjuntos difusos corresponde al “ indice de 
vulnerabilidad del puente”. 
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ABSTRACT 

The aim of this paper is to develop a simple model to assess vulnerability of bridges. This 
model was especially made for areas where current reports of seismic damages of bridges do 
not exist. The method is based on the identification of the most influential characteristics 
associated with the damage in a bridge under the action of an earthquake. The evaluation of 
these characteristics was carried out by means of the determination of nineteen parameters. 
A vulnerability degree and an importance value based on the opinion of experts was 
assigned to each of the 19 parameters. Due to the imprecision and subjectivity of the 
opinions, these data were processed by using fuzzy techniques. Once each one of the 
parameters was qualified, the qualifications of the parameters were related to their value of 
importance by means of a fuzzy weight average. The result of this relationship, according to 
techniques based on fuzzy sets, is the “ bridge vulnerability index". 



1 INTRODUCCI6N 



A1 ser los puentes una de las partes mas importantes dentro de la infraestructura 
viaria de un pais, surge la necesidad de estudios conducentes a la valoracion de su 
vulnerabilidad sfsmica. Ademas, no se debe ignorar que durante largos perfodos de servicio, 
los puentes existentes pueden estar expuestos a serios problemas sfsmicos y, en muchos 
casos, si el sismo es suficientemente fuerte, la estructura puede alcanzar incluso el colapso. 

Sin embargo, los estudios de vulnerabilidad son muy recientes, y debido a su 
importancia, se han convertido hoy en dfa en uno de los principales ffentes de investigation 
en ingenierfa sfsmica en todo el mundo. La mayorfa de los estudios de vulnerabilidad de 
puentes estan basados en datos obtenidos a partir de inventarios actualizados de las 
estructuras y de danos registrados ante efectos sfsmicos; estudios que solo son aplicables a 
las zonas de donde provienen los datos obtenidos. Por esta razon es muy importante 
proponer una metodologfa adaptable a regiones donde no se cuente con suficiente 
information sobre danos de puentes. 

El objetivo principal de esta investigation es presentar una metodologfa para evaluar 
la vulnerabilidad sfsmica de cada uno de los puentes de un area determinada. El metodo se 
basa en la definition de las caracterfsticas mas relevantes de un puente que pueden influir en 
su vulnerabilidad sfsmica y, basandose en la opinion de expertos, valorar la influencia de 
cada parametro. Las caracterfsticas citadas se han determinado mediante la election de 
diecinueve parametros, obtenidos de estudios realizados relativos al comportamiento 
sfsmico de puentes, experiencias posterremoto, estudios de modelos existentes y opiniones 
de expertos. Los parametros se clasificaron de acuerdo a cuatro condiciones de calidad, a las 
que se les asigno un grado de vulnerabilidad determinado y, a su vez, a cada parametro se le 
definio su valor de importancia. Estos grados de vulnerabilidad y valores de peso se 
identificaron a partir de opiniones de expertos. 

Dado que el modelo propuesto involucra opiniones de expertos y su posterior 
aplicacion se basa en la respuesta de un formulario donde se califica cada parametro, se 
decidio utilizar tecnicas basadas en conjuntos difusos. Dentro de los conjuntos difusos se 
utilizo el peso promedio ponderado difuso como medio para el calculo del indice de 
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vulnerabilidad sismica de puentes Este indice es la relation entre el grado de 
vulnerabilidad de cada parametro y su valor de importancia. 

El modelo se valido utilizando una muestra de trece (13) puentes, afectados por el 
sismo de Northridge del 17 de Enero de 1994, y que presentaban distintos grados de dano. 
Finalmente, la metodologia propuesta fue aplicada a algunos puentes de la provincia de 
Barcelona. 



2 PARAMETROS DEL MODELO 



La identification de los parametros mas influyentes en la evaluacion de la 
vulnerabilidad sismica de un puente, se efectuo en base a cuatro estudios realizados: 1) 
estudios sobre comportamiento sismico de puentes (Maldonado et al 1998 a y b), 2) 
metodologfas existentes de evaluacion de la vulnerabilidad sismica en puentes (Maldonado 
et al 1998 c), 3) experiences posterremoto y 4) opiniones de expertos. Con los puntos 1, 2 
y 3 se determinaron los parametros mas importantes en la determination de la 
vulnerabilidad sismica de puentes y con el punto 4, opiniones de expertos, se corroboro y 
complemento el trabajo realizado hasta el momento. 

A continuation, se expone cada uno de los cuatro estudios realizados . para la 
definition de los parametros que determinan la vulnerabilidad sismica de puentes en el 
modelo propuesto. 



2.1 Estudios sobre comportamiento sismico de puentes 



En los estudios de comportamiento sismico de puentes presentados en (Maldonado 
et al 1998 a y b), se estudian las caracterfsticas estructurales de los puentes que definen su 
respuesta ante la action sismica. Es asf como en este apartado tan solo se describen las 
conclusiones de los aspectos que son de interes en el presente trabajo. Para plasmar estas 
conclusiones, se estudio una pequena pero representativa muestra de puentes, la cual sirvio 
para definir algunos de los parametros que influyen en la respuesta dinamica de la estructura 
y que, por lo tanto, van a jugar un papel importante en el momento que se presente el dano 
en la estructura. 

La muestra de puentes de los que se sacaron algunos parametros, esta compuesta de 
seis puentes, cuatro corresponden a estructuras que en el momento de realizar su estudio 
estaban muy proximas a su construccion y dos estructuras ya existentes. Por ello, sobre las 
estructuras proyectadas pero sin construir se pudo realizar un analisis mas completo, de tal 
manera que las conclusiones sirvieran para las decisiones finales en su proyecto y 
construccion. En las estructuras existentes se evaluo el comportamiento de la estructura 
realizando analisis dinamicos que permitieron observar la influencia en la respuesta 
dinamica de sus elementos como, por ejemplo, las pilas, apoyos y juntas. 
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Los parametros con un mayor grado de influencia respecto a la respuesta sismica 
fueron los siguientes: existencia de articulaciones intemas, tipo de pila, tipo de cimentacion, 
sistema de apoyo de superestructura a subestructura, forma de la superestructura y 
regularidad longitudinal del puente. 



2.2 Metodologias existentes de evaluacion sismica de puentes 



De las diferentes metodologias existentes de evaluacion sismica que involucran 
modelos de la evaluacion de la vulnerabilidad sismica de puentes, se escogieron cinco 
(Kiremidjian y Basoz 1997, Ren y Gaus 1996, Dicleli y Bruneau 1996, Pezeshk et al 1993, 
ATC-6-2 1983), realizandose un estudio comparativo entre los distintos parametros 
considerados en sus modelo. 

Luego, se eligieron algunos parametros considerados en mas de un modelo, siendo 
estos: tipo de estribo, tipo de superestructura, tipo de subestructura, angulo de esviaje, tipo 
de cimentacion en pilas y estribos, historia de mejoramiento sismico del puente, ano de 
proyecto, armadura de las pilas, longitud de apoyos, existencia de articulaciones intemas, 
irregularidad en geometrfa y rigidez, condition de sitio, potencial de licuefaccion y tipo de 
apoyos. El valor de la aceleracion maxima en el sitio del puente no se selecciono dado que 
el modelo a proponer no pretende llegar hasta la evaluacion del dano de la estructura. 

2.3 Experiencias posterremoto 



De lo aprendido de sismos pasados, en lo referente a puentes, se pueden extraer las 
caracterfsticas que los hicieron vulnerables frente a la accion sismica. Para ello, se 
escogieron cuatro sismos los cuales causaron danos significativos en los puentes. Los 
sismos analizados son: el de San Fernando en 1971, Loma Prieta en 1989, Northridge en 
1994 y Kobe en 1995. 

Al comparar los danos causados en los puentes por cada uno de los cuatro sismos 
anteriores, se determinan las caracterfsticas (parametros) que hicieron que estas estructuras 
fuesen mas vulnerables frente a la accion del sismo: armadura en las pilas, especificaciones 
del proyecto, tipo de estribo, forma de la superestructura (angulo de esviaje), detalle de 
uniones, existencia de articulaciones intemas, longitud de apoyos y tipo de pilas. 

2.4 Opinion de expertos 



Basados en los tres puntos citados anteriormente se realizo una primera selection de 
diecisiete (17) parametros (ano de proyecto y construction del puente, tipo de 
superestructura, forma de la superestructura, existencia de articulaciones intemas, material 
de la superestructura, tipo de pila, tipo de cimentacion, material de las pilas, irregularidad 
longitudinal en geometrfa y rigidez, longitud de apoyo en pilas, tipo de suelo, tipo de estribo, 
longitud de apoyo en estribos, tipo de aparato de apoyo, estado de conservation del puente, 
procedimiento constructivo de la superestructura y procedimiento constructive de las pilas). 
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los cuales se dividieron, a su vez, en varias categonas; y, de esta manera, se pidio la opinion 
a cincuenta ( 50 ) expertos acerca de su consideracion, grado de vulnerabilidad y valor de 
importancia. De los cincuenta ( 50 ) expertos se recibieron veinticinco respuestas ( 25 ) y, con 
sus aportaciones, se redefinieron los parametros a considerar, quedando finalmente 
diecinueve ( 19 ) parametros. 

Los expertos en sus opiniones, consideraron que otros parametros podrfan ser 
importantes en la evaluacion de la vulnerabilidad sismica de puentes. Estos fueron 
analizados y reconsiderados en los ya propuestos, ya que muchos podi'an involucrarse en las 
categonas formuladas para cada parametro, considerandose algunos de ellos como nuevos 
parametros. 

Basados en los cuatro estudios realizados, anteriormente resumidos, se seleccionaron 
los siguientes parametros (Figura 1 ): 

(Ki) Ano de proyecto y construction del puente 

(K2) Tipo de superestructura 

(K 3 ) Forma de la superestructura 

(K 4 ) Existencia de articulaciones intemas 

(K5) Material de la superestructura 

{Ka) Tipo de pila 

(K7) Tipo de cimentacion 

(Kg) Material de las pilas 

(Kg) Irregularidad longitudinal en geometna o rigidez 

(K10) Longitud de apoyo en pilas 

(K11) Tipo de suelo 

(K12) Tipo de estribo 

(K13) Longitud de apoyo en estribos 

(K14) Tipo de aparato de apoyo 

(Kis) Estado de conservacion del puente 

{Ku) Procedimiento constmctivo de la superestructura (hormigon) 

(K17) Procedimiento constmctivo de las pilas (hormigon) 

(Kih) Potencial de licuefaccion 
( K19 ) Elementos no estmcturales 
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1. Ano de diseno y construccibn del puente 

2. Tipo de superestructura 

3. Forma de superestructura 

4. Existencia de articulaciones intemas 

5. Material de superestructura 

6. Tipo de pila 

7. Tipo de cimentacibn 

8. Material de las pitas 

9. Irreguiaridad longitudinal en geometria o rigidez 

10. Longitud de apoyo en pilas 



11. Tipo de suelo 

12. Tipo de estribo 

13. Longitud de apoyo en estribos 

14. Tipo de aparato de apoyo 

15. Estado de conservacibn del puente 

16. Procedimiento constructive de la 
superestructura (hormigon) 

17. Procedimiento constructive de las pilas (hormigon) 

18. Potencial de licuefaccibn 

19. Elementos no estructurales 



Figura 1 . Parametros del Modelo. 
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3 DETERMINACION DE LAS CONDICIONES DE C ALIDAD DE CADA 
PARAMETRO 

Siguiendo el mismo procedimiento realizado para la selection de estos parametros, 
de cada uno de ellos se consideraron tres o cuatro condiciones de calidad. Es asf como, 
basados en el estudio del comportamiento sfsmico de puentes, en las experiencias 
posterremoto y en la comparacion de metodologfas existentes de evaluacion sismica que 
involucran modelos de vulnerabilidad, se efectuo una primera subdivision de cada uno de 
los parametros identificados. Luego, cada parametro se dividio en tres o cuatro condiciones 
de calidad (A, B, C, D) (Maldonado et al. 2000 a, Maldonado 2000). Por ejemplo para el 
parametro que estudia la influencia del terreno sobre el cual se encuentra cimentado el 
puente, las divisiones impuestas fueron: A: Roca, B: Suelo granular denso o arcillas 
preconsolidadas ngidas, C: Arenas de densidad media (limosas), D: Arcillas medias a 
blandas. Para el parametro que evalua la influencia del tipo de superestructura las 
condiciones de calidad establecidas fueron: A: Puentes suspendidos, atirantados y puentes 
de un unico vano cuyo sistema corresponde a losa, vigas y sistemas de tablero en losa, vigas 
T, vigas cajon, vigas canal (artesa), vigas en celosfa, portico o arco, B: Puentes de dos o mas 
vanos con superestructura en sistema portico o arco, C: Puentes de dos o mas vanos con 
superestructura continua construida en sistema de losa, vigas y tablero en losa, vigas T, 
vigas cajon, vigas canal (artesa) o viga en armadura, D: Puentes de dos o mas vanos 
construidos en vigas simplemente apoyadas. 



4 REALIZACI6N DE ENCUESTAS A LOS EXPERTOS 



Una vez identificados los diferentes parametros que pueden influir en la 
vulnerabilidad de los puentes y los distintos aspectos a considerar dentro de cada parametro, 
se continuo con la determinacion de la calificacion de cada parametro. Para este fin se 
considero la opinion de expertos en el tema. 

Las opiniones fueron recopiladas mediante la realizacion de dos encuestas. La 
primera tenia como finalidad obtener las calificaciones de los parametros y, la segunda, la 
determinacion del valor de importancia de cada parametro dentro de la determinacion de la 
vulnerabilidad total del puente. Para las calificaciones de los parametros, cada uno de ellos 
se dividio en tres o cuatro condiciones de calidad y, de cada uno de ellas, se pregunto el 
grado de vulnerabilidad que se creia era el mas conveniente. 

Con la finalidad de obtener information fiable en las respuestas de las encuestas, se 
buscaron ihgenieros expertos en los temas de puentes y vulnerabilidad sismica. No siendo 
tarea facil, se decidio ampliar el ambito de expertos fuera de Espana, enviando las encuestas 
a otros paises de America como Colombia, Mexico, Ecuador y Venezuela. 

La election de una muestra de expertos fue una tarea laboriosa y delicada, ya que se 
pretendia que los especialistas contaran con un buen conocimiento del tema, una adecuada 
experiencia en el campo de trabajo y disponibilidad en el momento de responder las 
encuestas. La muestra total estuvo formada por cincuenta (50) expertos, el numero de 
encuestas enviadas a cada uno de los cinco paises fueron diferentes, dado que la base de 
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datos de los expcrtos con la que se contaba en el momento de enviar las encuestas no era 
homogenea en cantidad de expertos por pafs. 

Las encuestas fueron pensadas y elaboradas buscando que fueran breves, concisas y 
claras. Se construyeron en forma de tablas para facilitar la forma de respuesta y, al mismo 
tiempo, para organizar la informacion de tal manera que el encuestado no perdiera el sentido 
de la pregunta. 

Las encuestas fueron, como se ha citado, dos: la primera tenia como finalidad el 
comparar el grado de vulnerabilidad (“ nada ”, “ poco ”, “medianamente", “ muy ” y 
“ absolutamente ” vulnerable) de cada una de las clasificaciones (A, B, C y D) para un 
determinado parametro y asignar a cada division un grado de vulnerabilidad; la segunda, 
tenia como proposito el conocer la importancia relativa que existe entre los diferentes 
parametros. En la Figura 2 y 3 se presenta el modelo de la primera y segunda encuesta para 
el primer parametro. 



PRIMERA ENCUESTA 

PROPOSITO DEL CUESTIONARIO: El proposito de este cuestionario ej 
comparar el grado de vulnerabilidad de cada una de las clasificaciones A, B,C y f 1 
dentro del parametro K/ respectivo. 

INSTRUCCIONES DE LLENADO: A continuacion se presenta una tabla la cua 
contiene los parametros K, que se cree influyen en la determinacion de le 
vulnerabilidad sismica de un puente, y -sus correspondientes clasificaciones A, B 
C y D. Cada una de estas clasificaciones debera ser llenada de acuerdo a l< 
siguiente escala de valores dado el grado de vulnerabilidad que considere para cad£ 
clasificacion: 



NADA Vulnerable 


0 


POCO Vulnerable 


3 


MEDIANAMENTE Vulnerable 


5 


MUY Vulnerable 


7 


ABSOLUTAMENTE 


9 


Vulnerable 





1 Ki 


- ANO DE PROYECTO Y CONSTRUCCION DEL PUENTE 




KI 


Despues de 1995. 




1 




Entre 1975 y 1994. 




1 




Entre 1968 y 1974. 




1 




Antes de 1968. 




u 



Figura 2. Modelo de la primera encuesta para el parametro Kj. 
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SEGUNDA ENCUESTA 

PROPOSITO DEL CUESTIONARIO: La finalidad de este cuestionario es 
conocer la importancia relativa que existen entre los parametros K,_ 

INSTRUCCIONES DE LLENADO: A continuation se presenta una tabla con los 
diecisiete parametros solamente. Se debera llenar cada fila con el valor que cree 
conveniente de acuerdo a la importancia que considere tiene el parametro er 
comparacion con los restantes. La escala de calificacion es de 0 a 10, pudiendoss 
llegar a parametros que considere igualmente importantes y por lo tanto tendran U 
misma calificacion. 



NADA Importante 


0 


MAXIMA Importancia 


10 



■ 


PARAMETRO 






Ano de Proyecto y construction 


i 



Figura 3. Modelo de la segunda encuesta para el parametro Kj. 



La buena aceptacion de las encuestas por parte de los expertos se reflejo en el 
numero de las recibidas y en las valiosas aportaciones que muchos de los especialistas 
anexaron como observaciones. 



Durante algunos meses se estuvo a la espera de la respuesta y una vez recibidas las 
mfnimas deseadas, se dio paso al procesado de los datos recibidos. Los datos una vez 
adquiridos, se tomaron en su totalidad como aceptables, dada la fiabilidad de las respuestas 
por parte de los expertos. De esta manera, se utilizaron todas las respuestas para construir las 
funciones de pertenencia necesarias para el desarrollo del modelo. 



De los cincuenta expertos a los que se les envio las encuestas, correspondientes a 
cinco paises, se obtuvo una buena cantidad de respuestas. Veinticinco expertos (50%) 
enviaron sus opiniones. El visto bueno de las encuestas fue aceptable en todos los paises, 
pero dado el mayor numero de encuestas enviadas a Colombia, fue de este pais de donde se 
obtuvo la mayor cantidad de respuestas (Figura 4). 




a. Distribution de la muestra total de expertos a los que se enviaron encuestas 
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Venezuela 
Mexioo 8% 




Ecuador 



Colonnbia 

52 % 



12% 



b. Distribution de muestra de opiniones de expertos recibidas 
Figura 4. Muestra de expertos. 



5 TRATAMIENTO DE LA INFORMACION FOR CONJUNTOS DIFUSOS 



Dado que no se contaba con informacion de danos en puenteS ocurridos en sismos 
pasados, fue necesario recurrir a la opinion de expertos para obtener los datos necesarios 
para definir los grados de vulnerabilidad y los valores de importancia. Luego, al provenir la 
informacion de opiniones esta contiene subjetividad, imprecisiones e incertidumbres; razon 
por la cual una buena altemativa de manejar la informacion es mediante la teoria de 
conjuntos difusos (Maldonado et al. 2000 b). 



5.1 Que son los conjuntos difusos? 



El concepto de conjunto difuso lo introduce Zadeh en 1965 (Zadeh 1965), como un 
intento de veneer la severidad de la teoria clasica de conjuntos y poder reunir proposiciones 
que, por la naturaleza de lo que representan contienen incertidumbre, imprecision, 
ambigiiedad, generalidad, errores, aproximaciones, borrosidad. La teoria difusa es una teoria 
matematica, que se llama difusa porque recoge un aspecto de incertidumbre. Lo difuso es la 
ambigiiedad que puede ser encontrada en la definicion de un concepto o en el significado de 
una palabra. Por ejemplo, la incertidumbre en las expresiones como: “puente levemente 
dafiado ”, “grieta considerable", “estructura funcionalmente obsoleta “, “alt a seguridad' o, 
“estructura poco vulnerable”. 

Hablando en terminos mas simples, la “Teoria de Conjuntos Difusos” parte de la 
teoria clasica de conjuntos, anadiendo una funcion de pertenencia al conjunto, p, 
generalmente definida como un numero natural entre 0 y 1. Asl, se introduce el concepto de 
conjunto difuso asociado a un determinado valor lingiifstico, definido por una palabra, 
adjetivo o etiqueta lingiustica A. Para cada conjunto difuso se define una funcion de 
pertenencia o inclusion Pa(x), que representa el grado en que la variable x esta incluida en el 
concepto representado por la etiqueta A. 

Los conjuntos difusos permiten reunir objetos o sucesos por el valor de una cierta 
magnitud. Por ejemplo, las edificaciones para vivienda pueden ser reunidas por su edad. 
Asi, si definimos el conjunto clasico de edificaciones viejas como las que tienen mas de 30 
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anos, resulta que una de 29 anos no serfa vieja, pero otra de 31 si que lo serfa, siendo la 
diferencia entre ellas de tan solo dos anos. Por lo tanto, una descripcion del conjunto de las 
edificaciones viejas en terminos de conjuntos difusos parece mas adecuada en este tipo de 
casos. 



La teorfa de los conjuntos difusos trabaja con la cuantificacion de los significados de 
las palabras en graficas, dentro del marco de trabajo de la teorfa de conjuntos. Es un intento 
de expresar “alto” o “ viejo ” por medio del concepto de conjuntos. Desde un punto de vista 
practico, el concepto de conjunto no es necesario para la cuantificacion de la ambigiiedad, 
pero es posible aumentar el rango de utilidad trabajando dentro del marco de trabajo de la 
teorfa de conjuntos. Esto se debe a que la teorfa de conjuntos es un concepto muy basico y 
esta relacionada con todos los campos de las matematicas contemporaneas. 



5.2 Representation de los conjuntos difusos 

Si X es una coleccion de objetos, denotados generalmente por x, X={xj, X 2 x n f, 

luego un subconjunto difuso A en X es un conjunto de pares ordenados: 

A = {p A (x)\x, x e X} (1) 

Donde p A ( x) es llamada la funcion de pertenencia o grado de pertenencia (tambien 
grado de compatibilidad o grado de verdad) de x en A, es decir Pa(x) define el grado en el 
cual el elemento jc del conjunto X esta incluido en el subconjunto A. Es asf como la funcion 
de pertenencia es la que define el conjunto difuso. 

El rango de la funcion de pertenencia puede ser un subconjunto de numeros reales 
no negativos, aunque sea de practica general que la funcion de pertenencia este definida 
entre 0 y 1 como Pa(x) ■' X— >[1,0]. 

El grado de inclusion, algunas veces, se denomina la magnitud ( extend ) o el grado 
(grade). Por ejernplo, el grado de pertenencia del elemento x en el subconjunto A se expresa 
mediante: 



p A (xf) = \, /i A (x 2 ) = 0.8, p A (x 3 ) = 03, fi A (x 4 ) = 0 

Donde p es la funcion de pertenencia, proporcionando el grado de pertenencia, un valor 
que varfa de 0 hasta 1. El subindice de p, A, expresa que p,\ es la funcion de pertenencia de 
A. 



5.3 Definicion de funcion de pertenencia 

La teorfa de los conjuntos difusos define el grado en el cual el elemento x del 
conjunto X esta incluido en el subconjunto A mediante la llamada funcion de pertenencia 
( membership function), Pa(x). Dicha funcion es un conjunto de numeros ordenados si la 
variable es discreta, o una funcion continua si no lo es. Como se ha mencionado 
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anteriormente, repetida veces, el valor Pa(x) indica el grado en que el valor x de la variable 
X esta incluida en el concepto representado por la etiqueta A. 



5.4 Concepto de variable lingiifstica 

Se denomina variable lingiifstica a aquella que puede tomar por valor terrrunos del 
lenguaje natural, como edad, peso, altura, dano, capacidad, vulnerabilidad.... nada, poco, 
mucho, positive, negativo, bueno, malo, etc., las cuales a su vez son las que hacen el papel 
de etiquetas en un conjunto difuso. 

Aun cuando la finalidad primordial de este concepto es expresar de manera formal el 
hecho de que pueda asignarse como valor de una variable palabras tomadas del lenguaje 
natural, no obstante, a una variable lingiifstica podran asignarse valores numericos. Asf, en 
una expresion como el dano estructural es severo, la variable dano estructural debe ser 
entendida como una variable lingiifstica, pues se le asigna como valor el conjunto difuso 
severo, pero, ademas, esta variable puede tambien tomar valores numericos como el indice 
de dano es mayor de 1. 



5.5 Variables lingiusticas utilizadas 

En el algoritmo de desarrollo para el calculo del fndice de vulnerabilidad mediante la 
utilizacion de opiniones de expertos fue necesario definir adecuadamente las variables 
lingiifsticas para describir los diferentes grados de vulnerabilidad. 

Investigaciones sobre evaluacion del dano de estructuras o evaluacion del estado de 
estructuras afectadas por el sismo, en las cuales el uso de la forma lingiifstica se maneja 
mediante teorfa difusa, han servido de base para la definition de las variables utilizadas en 
este estudio. Estas variables lingiifsticas han sido utilizadas para calificar cada una de las 
condiciones de calidad impuestas a cada parametro. A su vez, se han empleado tambien para 
asociarle a cada funcion de pertenencia del fndice de vulnerabilidad calculado un cierto 
grado de vulnerabilidad. 

Asf pues, las variables lingiifsticas naturales seleccionadas para calificar los grados 
de vulnerabilidad son “ nada ”, “ poco ”, “ medianamente ”, “muy” y “absolutamente” 
vulnerable. Cada variable requirio la asignacion de un valor numerico, escogiendose 
arbitrariamente una escala comprendida entre 0 y 9 para identificar cada grado de 
vulnerabilidad. 



Los valores numericos asignados a cada variable para la materialization de las 
encuestas fueron: 

0 NADA Vulnerable 

3 POCO Vulnerable 

5 MEDIANAMENTE Vulnerable 
7 MUY Vulnerable 

9 ABSOLUTAMENTE Vulnerable 
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Posteriormente, los valores citados se representaron por medio de conjuntos 
difusos compuestos de valores numericos. En esta investigation, los conjuntos difusos que 
representan las variables relativas al grado de vulnerabilidad, se asumen como: 



Nada 

Poco 

Medianamente 

Muy 

Absolutamente 



{ 1 | 0 , 0 | 0 } 

{ 0 | 0 , 1 | 3 , 0 | 5 } 
{0|3 , 1|5 ,0|7} 
{0|5 , 1|7 , 0|9} 
{0|7 , 1|9 } 



Siendo de la forma: 



y = {fMx)\x} 



( 2 ) 



donde /j(x) denota el grado de pertenencia de x, definiendo jc como el universo de la 
vulnerabilidad del puente en los conjuntos difusos. 



Para la definition de las funciones de pertenencia de las variables que definen el 
grado de vulnerabilidad del puente se toma solamente el universo formado por los valores 0, 
3, 5, 7 y 9. El involucrar mas elementos en el universo ayudana a representar mejor la 
cantidad de interes, no obstante, tambien involucrarfa una mayor manipulation aritmetica en 
los calculos. Por ello, se eligio tomar el conjunto difuso con solo los cinco elementos 
citados. 



La representation grafica de las funciones de pertenencia de las variables linginsticas 
utilizadas en este estudio es la siguiente: 




Figura 5. Funciones de pertenencia de las variables linginsticas “nada", “poco”, 
“medianamente” , “muy” y “absolutamente” vulnerable. 

5.6 Construccion de las funciones de pertenencia 

Para la construccion de las funciones de pertenencia se aprovecho la information de 
las respuestas de las encuestas obtenidas de 25 expertos de cinco pafses: Espana, Colombia, 
Ecuador, Venezuela y Mexico. Los valores de pertenencia se calcularon en base al numero 
de respuestas favorables de cada clasificacion particular (Maldonado et al. 2000 a. Tee et 
al. 1988). 
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Dada la subjetividad y ambigiiedad de las opiniones de los expertos, muchas de las 
respuestas estuvieron encasilladas en los grados de vulnerabilidad de “poco" a 
“ medianamente ” vulnerable. De ahf, que las funciones de pertenencia construidas 
inicialmente se corrigieron de acuerdo con la forma de calculo de un fndice de 
vulnerabilidad sfsmica expuesto posteriormente. 

Corregidas las funciones de pertenencia de las calificaciones, estas se incorporaron 
dentro del algoritmo para el calculo del fndice de vulnerabilidad sfsmica de puentes. Como 
ejemplo, a continuation se presentan las funciones de pertenencia para las clasificaciones de 
uno de' los diecinueve parametros, el Kz denominado “ tipo de superestructura'. 




Puentes suspendidos, atirantados y puentes de un unico vano. 

Puentes de dos o mas vanos con superestructura en portico o arco. 
Puentes do dos o mas vanos con superestructura continua. 

Puentes de dos o mas vanos construidos en vigas simplemente apoyadas. 



Figura 6. Funcion de pertenencia del grado de vulnerabilidad del parametro “tipo de 

superestructura ” 




Figura 7. Funcion de pertenencia del valor de importancia para el parametro “ Tipo de 

superestructura” > 

De esta manera construidas y corregidas las funciones de pertenencia se llego a la 
definition de los grados de vulnerabilidad de las condiciones de calidad de los parametros 
que influyen en la vulnerabilidad sfsmica de los puentes y de sus valores de importancia, a 
partir de las opiniones de expertos y mediante la utilization de “conjuntos difusos 
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5.7 Numeros difusos e intervalos difusos 

La idea basica de un numero dijuso es que debe representar algun valor numerico, el 
cual ha sido definido imprecisamente. Sin embargo, siendo mas especificos, se debe definir 
a los numeros e intervalos difusos como: 

Un numero difuso A es una cantidad difusa p A \ R — » [0,1] que satisface las 
siguientes condiciones : 

a) A debe ser convexo. Una idea basica conectada con que un numero difuso debe 
representar algun valor numerico determinado imprecisamente es que la funcion de 
pertenencia de un numero difuso no debe tener distintos puntos maximos, es decir que los 
numeros difusos deben ser “ convexos \ 

b) Pa esta normalizado. Existe un x que pertenece a los R con p A (x)=l (Este punto x se 
denomina “ valor medio ” de A) 

Esta definicion implica que los numeros difusos son simplemente conjuntos difusos 
reales A, para los cuales los a-cortes A a -{xeX ;p A (x)>a} son intervalos ceirados y acotados 
de la recta real. 



5.8 Operaciones aritmeticas con numeros e intervalos difusos 

Ademas de la aritmetica conectiva y de la aritmetica ponderada, existe un tercer tipo 
de aritmetica ligada a las cantidades o a los numeros boirosos. Para calcular con numeros e 
intervalos difusos, primero se necesita la definicion de las operaciones aritmeticas basicas 
para ellos. Para ello, se aplica el principio de extension (Zadeh 1965). Luego la suma, 
multiplicacion y la division de dos numeros o intervalos difusos se determina como: 

A partir de la definicion de los conjuntos Ay B como: 

A = {ftfi) | i} (3) 

B = ffMj)\jJ ( 4 ) 

donde fi(i) y jj(j) son los grados de pertenencia de los elementos i y j, respectivamente. 

Las operaciones citadas se expresan como: 



a) Adicion 



A + B = max{min(n A (i),V B (j))\[i + j ] } 



( 5 ) 



b) Producto 



c) Division 



A B = max{min(n A {i),jx a (y))|[i • ;] } 



( 6 ) 
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A® B - max{min(fx A (i),fi B (;'))|[i © j]} 



(7) 



Los conjuntos difusos resultantes de estas definiciones, son siempre intervalos 
difusos si A y B son intervalos, y son numeros difusos si A y B son numeros difusos. 

Estas operaciones aritmeticas para numeros e intervalos difusos, generalizan las 
operaciones aritmeticas usuales para los numeros reales y tambien las operaciones de los 
intervalos aritmeticos. Muchas de las leyes que se cumplen para la aritmetica de numeros 
reales, se generalizan a estas operaciones para intervalos difusos, aunque no todas ellas. 



5.9 Distancia entre conjuntos difusos 

Una vez realizados calculos que incorporan operaciones como la suma, 
multiplication y la division difusas se obtiene un conjunto difuso. De tal manera que para 
interpretar mejor lo obtenido, generalmente, se requiere transformar el conjunto difuso 
resultante a una expresion del lenguaje natural. Dicha transformation puede ser un proceso 
relativamente simple si se involucra la determination de la distancia entre un conjunto 
difuso dado (resultante) y los conjuntos que representan cada una de las posibles 
expresiones del lenguaje natural (variables lingiifsticas). 

En general, la distancia entre dos conjuntos difusos A y C se puede definir utilizando 
diversas medidas, siendo las mas frecuentes las siguientes: 

a) Hamming 

n 

Distancia ( A, C) = ^ (fi A (x) - fi c (x)) (8) 



b) Euclfdea 



Distancia (A,C) = 



n 



y. (■*) Be (•*)) 




(9) 



c) Minkowski 



Distancia ( A, C) = 



X 



w 



(H A (x) - n c (x)) w 



con we[l,«] 



( 10 ) 



5.10 El problema de la desfusificacion 

Los sistemas basados en sistemas difusos, generalmente, originan resultados 
apoyados en sus reglas lingiifsticas, que corresponden a conjuntos difusos. Pero en muchas 
de las aplicaciones, los resultados que se requieren son precisamente los “no difusos”. Para 
transformar los resultados difusos en no difusos, se ha desarrollado la desfusificacion. 
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Siendo el desfusificador la funcion que transforma un conjunto difuso, el cual normalmente 
corresponde a la salida de un dispositivo de inferencia difusa, en un valor no difuso. Para 
esta tarea existen diversos metodos, que se pueden resumir en los procedimientos siguientes: 

a) Procedimiento maximo. 

La forma mas simple es el llamado procedimiento-maximo. En el, se supone que la 
funcion de pertenencia pg de la salida difusa B tiene un punto maximo simple, tomando a 
continuation el valor de desfusifcacion (yo(B)) el argumento (abscisa) en este punto 
maximo: 



y 0 (B) = &rg max {p B (y)\y eY } (11) 

No obstante, esta presuncion es poco especifica. Si Pb tiene varios valores maximos, 
se tendra mas de un punto maximo. En estos casos, se puede empezar el proceso con el 
conjunto B inax de todos los argumentos de puntos maximos de Pb 

B max = { yeY | p(y) = max p B (z) } (12) 

La idea es que Bmax es el conjunto de decisiones optimas. Debiendose, finalmente, 
escoger uno de ellos. 

Una de las formas de escoger una decision final es el seleccionar un elemento de 
Umax aleatoriamente, de tal forma que se presupone que todos los elementos de B nuLX 
corresponden a decisiones igualmente buenas. La interpretation de que los puntos que estan 
fuera de B nmx no son optimos se entiende, frecuentemente, como una indication de que yo(B) 
no debe corresponder a la region del borde de B tnax , preferiblemente de la mitad de ella. 

Si se supone, adicionalmente, que el universo del discurso Y de la variable de salida 
v, es un conjunto de numeros, Y qR, otra propuesta es tomar el valor medio de B nulx . En el 
caso de que Bmax sea un conjunto finito, el valor medio vendra dado por: 

y o( B ) = ~^ < 13 > 

siendo N , el numero de elementos de B max (cardinal). Pero se debe tener cuidado aquf, ya 
que un valor medio de B nmx puede ser una decision preferible en el caso que B mix sea un 
intervalo en Y, pero siempre se debe reflexionar sobre la decision finalmente tomada. 

b) Centro de gravedad. 

Otra desventaja de cualquier tipo de procedimiento maximo es que se ignora toda la 
informacion concemiente a la funcion de pertenencia ps fuera de Bmax. Una forma de 
considerar toda la funcion de pertenencia es utilizar el procedimiento del centro de 
gravedad (COG), el cual corresponde a una estrategia de valor medio. Este valor, que 
determinayoffl), se toma como la media de los pesos de todo el supp(B) (los argumentos del 
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conjunto difuso cuyas funciones de pertenencia son mayores que cero) como: 



card(B) ^ 



(14) 



siendo card (B) la suma de las funciones de pertenencia ^(y)- 



6 evaluaci6n del Indice de vulnerabilidad sismica difuso 



Una vez procesada la information como conjuntos difusos las aritmeticas difusas 
fueron empleadas para relacionar las calificaciones de cada parametro y sus respectivos 
valores de importancia, por medio del “peso promedio”. Y de esta manera esta expresion 
fue utilizada como un “ indice de vulnerabilidad sismica de puentes”, con el cual se pueda 
evaluar la vulnerabilidad sismica de puentes y; el cual se expresa como (Maldonado et al. 
2000 a, b, Maldonado 2000): 



£ W'K, 

I -V. 1=^9 (IS) 

i=l 



Donde L V.j es el indice de vulnerabilidad sismica del puente j, Ki es una medida del 
grado de vulnerabilidad de la categorfa del parametro i. Los valores Wi son una medida de la 
opinion de la importancia asociada al parametro i respecto a los demas parametros. 



7 ALGORITMO PARA CALCULAR EL INDICE DE VULNERABILIDAD 



El calculo del peso promedio difuso para la evaluacion del fndice de vulnerabilidad 
sismica de puentes se realizo de la siguiente forma: 

1) Traslacion de las variables de calificacion de cada una de las condiciones de calidad de 
los diecinueve parametros a conjuntos difusos, mediante la definicion de las 
correspondientes funciones de pertenencia. Estas funciones se construyeron a partir de las 
opiniones de expertos. 

2) Conversion de las variables de importancia de cada parametro a conjuntos difusos 
mediante las funciones de pertenencia. De igual manera que en el punto anterior, estas 
funciones de pertenencia se elaboraron en base a la opinion de expertos. 

3) Combinacion de las variables de calificacion y las variables de importancia difusas 
para obtener un conjunto difuso que represente el sistema completo, mediante la expresion 
75. Esta combinacion da como resultado un conjunto difuso llamado “ indice de 
vulnerabilidad ’. 
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4) Transformacion del conjunto difuso resultante, “indice de vulnerabilidad”, a una 
expresion lingiifstica como u nada”, “poco”, “medianamente ”, “ muy ” y “ absolutamente ” 
vulnerable. La transformacion de un conjunto difuso a una expresion lingiifstica natural es 
un proceso relativamente simple, el cual involucre la determinacion de la distancia del 
conjunto difuso resultante a cada uno de los conjuntos difusos que representan las 
expresiones lingiiisticas citadas. La minima distancia entre el conjunto difuso resultante y 
los conjuntos difusos que representan las variables lingiiisticas, es la expresion lingiifstica a 
la que se le asocia el fndice de vulnerabilidad calculado. 

5) Asociacion a la variable lingiifstica del fndice de vulnerabilidad calculado de un valor 
numerico no difuso. Esto se ha hecho mediante el proceso de desfusificacion. 



8 FORMA DE APLICACI^N DEL METODO 



Para la aplicacion del modelo, es necesario contar con el inventario previo de la 
muestra de puentes. Este inventario se realiza mediante el “formulario de levantamiento de 
la vulnerabilidad sismica de puentes El formulario contiene todos los parametros 
necesarios para la determinacion del mdice de vulnerabilidad del puente respectivo. Tan 
solo se requiere que el evaluador asigne a cada parametro su respectiva calificacipn. Una 
vez obtenidas todas las calificaciones (A, B, C y D) de cada uno de los diecinueve 
parametros analizados del puente, se procede al calculo del “mdice de vulnerabilidad”, I.V., 

El formulario de levantamiento de la vulnerabilidad sismica constituye una 
description completa de los datos necesarios para obtener el “indice de vulnerabilidad \ En 
base a la experiencia de los expertos y al estudio de los autores del metodo, como se expuso 
anteriormente, se han identificado diecinueve parametros, considerados los mas relevantes e 
influyentes en el dano que sufrirfa un puente sujeto a la action sismica. Algunos de estos 
parametros son puramente subjetivos, mientras que los restantes se obtienen a partir de 
mediciones o calculos simplificados que aseguran versatilidad y agilidad a la evaluation del 
fndice de vulnerabilidad. Cada parametro se clasifica en el formulario mediante una escala 
de A a D en forma decreciente, conforme la calidad del parametro vaya tambien 
decreciendo. 

Para la construction del formulario, se tomo como base algunos ya existentes para 
inventarios de puentes y, a partir de ellos y de los datos necesarios para el calculo del fndice 
de vulnerabilidad, se construyo el formulario de levantamiento sfsmico de puentes. El 
formulario se desarrolla en dos partes. En la primera, se detalla la information general de los 
puentes, tal como su denominacion, la localizacion geografica, la carretera soportada por el 
puente, la carretera que cruza bajo el puente, el propietario, el tipo de paso, los ensanches, 
las ampliaciones y lo mas importante, la calificacion de los parametros que influyen en su 
vulnerabilidad. La segunda parte del formulario corresponde a los “comentarios” , y se 
elaboro con la finalidad de servir de gufa en la respuesta a las preguntas del formulario. 
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Se intenta que el formulario sea facil de cumplimentar, partiendo de datos de 
inventarios existentes o de information que pueda adquirirse facilmente a partir de la 
inspeccion visual del puente. En la Figura 8 se presenta parte del formulario de 
levantamiento de la vulnerabilidad sfsmica de puentes. 






1. LOCAL! ZACION 

Denominacion 

P.K. 

Coordenadas x _ 



Termioo Municipal _ 

Provincia : 

Poblacion anterior _ 

Pobiacion posterior 
Propietario : 






1. ANO DE PROYECTO DE CONSTRUCCION DEL PUENTE 


A 








[ ] 




B 


Entre 


1 Q7R u 1 QQ4 




[ 3 




.....C 


Entre 


1 968 v 1 974 




[ ] 




D 


Antes 


de 1968 




[ ] 



2. TIPO DE SUPERESTRUCTURA 


A Puentes 


suspendidos, puente atirantados, puentes de un unico vano 


1 1 

B Puentes 


de dos o mas vanos en sistema portico o arco 


1 1 

C Puentes 


de dos o mas vanos con superestructura continua 


1 1 

D Puentes 


de dos o mas vanos en vigas simplemente apoyadas 


1 1 
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3. FORMA DE LA SUPERESTRUCTURA 




rectos (cero grados de esviaje) 


[ 1 

B Puentes 


esviados menos de veinte grados 


1 1 


esviados entre veinte y cuarenta y cinco grados 


[ ] 


con angulo de esviaje mayor a cuarenta y cinco grados. 




Puentes curvos 


[ ] 



Figura 8. Parte del formulario de levantamiento de la vulnerabilidad sfsmica de puentes. 



9 PROGRAMA INDICE PARA LA EVALUACION DEL INDICE DE 
VULNERABILIDAD SISMICA DE PUENTES 



Construidas las funciones de pertenencia para cada una de las tres o cuatro 
clasificaciones de cada uno de los parametros y determinadas las funciones de pertenencia 
de los valores de importancia, se calculo el “indice de vulnerabilidad" utilizando la 
expresion de peso promedio difuso de la ecuacion 15. 

El calculo de este fndice involucra una serie de operaciones aritmeticas con 
conjuntos difusos cofno la suma, el producto y la division. A su vez, en la transformation de 
los conjuntos difusos a una variable linginstica se incorpora el calculo de la distancia 
minima entre conjuntos difusos, todo ello sin obviar el control que se debe realizar a cada 
operacion aritmetica para que cumplan con la propiedad de convexidad y normalization. El 
programa INDICE se elaboro con el fin de incorporar todas las operaciones necesarias para 
el calculo del peso promedio difuso definiendo el fndice de vulnerabilidad. Las funciones de 
pertenencia, ya construidas, se convierten en la base de datos del programa y con la 
respuesta al formulario de levantamiento del puente, se escogen las funciones respectivas y 
se realiza el calculo del fndice. 

El programa INDICE se estructura en tres modulos basicos y una base de datos 
(Figura 9). Los tres modulos son: 1) Initialization del analisis, 2) Evaluacion del fndice de 
vulnerabilidad y 3) Analisis de los resultados. Estos modulos se alimentan a partir de la 
base de datos que contiene los datos de los puentes y las funciones de pertenencia requeridas 
para el calculo del fndice de vulnerabilidad. A continuation, se expone brevemente el 
algoritmo desarrollado a partir de estos tres modulos: 

1 ) Inicializa cion del analisis. En este modulo se define la muestra de puentes a trabajar. 

2) Evaluacion del fndice de vulnerabilidad para el puente P(j). El calculo del fndice 
de vulnerabilidad esta incorporado en este modulo. Aquf se lee el formulario de 
levantamiento del puente P(j), se hacen las respectivas operaciones aritmeticas difusas 
necesarias para calcular el fndice de vulnerabilidad. Efectuada cada operacion, el programa 
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INDICE controla que se cumplan las propiedades de convexidad y normalization. En caso 
de no cumplirse, corrige la funcion antes de continuar con los siguientes calculos. 

Este modulo esta alimentado a partir de la base de datos que contiene los resultados del 
formulario de levantamiento de cada puente de la muestra, las funciones de pertenencia de 
cada una de las calificaciones de los diecinueve parametros y las funciones de pertenencia 
de sus pesos. 

En el momento que INDICE lee los datos del puente P(j), selecciona las funciones de 
pertenencia respectivas a las calificaciones asignadas a cada parametro. 

3) Analisis de resultados del puente P(j). Una vez se ha calculado el mdice de 
vulnerabilidad sismica del puente P(j) mediante la ecuacion 1, este corresponde a un 
conjunto difuso y se expresa mediante una funcion de pertenencia. Razon por la cual es 
interesante y necesario interpretar el resultado obtenido para que asi sea de facil manejo y 
uso posterior. Este proceso dentro de la teorfa de los conjuntos difusos se llama 
desfusificacidn. El programa INDICE tiene dos opciones para realizar este proceso, una es 
transformar el conjunto difuso en una variable lingiiistica natural como “ nada’\ “poco”, 
“medianamente ” , “ muy ” y “ absolutamente ” vulnerable. Esto se efectua por medio de la 
minima distancia entre el conjunto difuso de mdice de vulnerabilidad y los conjuntos 
difusos de grados de vulnerabilidad. La segunda option del programa INDICE , es asignarle 
un valor numerico al conjunto difuso calculado. Esto lo hace el programa escogiendo el 
valor maximo de la funcion de pertenencia, es decir el valor de abscisa que le corresponde 
una ordenada de valor 1 . 

Este modulo requiere de la base de datos, la cual contiene las funciones de pertenencia de 
los diferentes grados de vulnerabilidad (“nada”, “ poco ”, “ medianamente’\ “muy” y 
“ absolutamente ” vulnerable). 



Base de Datos. En lo referente a la base de datos, el programa INDICE necesita de las 
respuesta de los formularios de levantamiento de cada uno de los puentes de la muestra. A 
su vez, el programa INDICE requiere leer las funciones de pertenencia tanto de las 
calificaciones de los parametros como de sus valores de importancia. Las funciones de 
pertenencia de los pesos o valores de importancia permanecen constantes a lo largo del 
calculo del mdice de vulnerabilidad de varios puentes. Las funciones de pertenencia de las 
calificaciones de los parametros vanan de acuerdo a la respuesta del formulario del puente 
en cuestion. 
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Figura 9. Algoritmo del programa INDICE . 
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10 CONCLUSIONES 

En este trabajo se ha expuesto un modelo para calcular la vulnerabilidad sfsmica de 

puentes a gran escala y se ha presentado el desarrollo de la metodologfa propuesta en todas 

sus facetas. Las principales observaciones y conclusiones son: 

1. El metodo propuesto constituye una description completa de la information necesaria 
para obtener un cuantitativo de la vulnerabilidad sfsmica de puentes llamado “ indice de 
vulnerabilidad'. El calculo del fndice se basa en la identification de diecinueve 
parametros, considerados los mas relevantes e influyentes en el dano que sufrira un 
puente bajo la action del sismo. Definiendo el “ mdice de vulnerabilidad" como la 
relation existente entre el grado de vulnerabilidad de cada parametro y su respectivo 
valor de importancia. 

2. Debido a que el modelo considera la opinion de expertos, para definir tanto el grado de 
vulnerabilidad de los diecinueve parametros como el valor de importancia de cada 
parametro, dentro del calculo de la vulnerabilidad total del puente, se decidio trabajar con 
conjuntos difusos, los cuales consideran y manejan de la mejor forma posible la 
subjetividad e imprecision de las respuestas de los especialistas. De esta manera, el 
“indice de vulnerabilidad' propuesto, se calcula a traves de un peso promedio difuso, el 
cual relaciona la medida del grado de.vulnerabilidad de cada parametro y sus respectivos 
valores de importancia. A1 provenir los resultados de las encuestas de opiniones de 
expertos, era de esperar que estos resultados fueran subjetivos y con incertidumbres y 
ambiguedades. 

3. El modelo busca ser de un ambito de aplicacion general; sin embargo, su utilizacion en 
Espana y America deben tener en cuenta condiciones de calidad particulares impuestas 
en unos pocos parametros. Esto se ideo desde el principio y, por ello, al categorizar cada 
parametro, se busco que fuera factible su utilizacion en cualquier parte. Sin embargo, el 
primer parametro definido “ano de especificacion de proyecto y construccion" del puente 
varfa de acuerdo a la normativa del pafs donde corresponda dicho puente. Por esto, este 
fue el unico parametro al que se le aplicaron dos tipos de clasificaciones, una para Espana 
y otra para los pafses de America en general. 

4. Los parametros considerados imprescindibles en el momento de la evaluation de la 
vulnerabilidad sfs mi ca del puente, segun las funciones de pertenencia, son el Ki (ano de 
proyecto y construction del puente), K 3 (forma de la superestructura), K 4 (existencia de 
articulaciones intemas), K 9 (irregularidad longitudinal en geometrfa o rigidez), Kjo 
(longitud de apoyo en pilas), Ku (tipo de suelo), K 13 (longitud de apoyos en estribos) y 
Kjs (estado de conservation del puente). 

5. Con el modelo propuesto, se obtiene un calificativo lingufstico y numerico del fndice de 
vulnerabilidad de un puente ademas de contarse con la funcion de pertenencia de dicho 
fndice. El modelo califica al puente entre nada, poco, medianamente, muy y 
absolutamente vulnerable. A su vez, le asigna un valor de 0 a 9 del fndice de 
vulnerabilidad, ademas de la funcion de pertenencia. 
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6. En la mayorfa de los casos los datos necesarios para aplicar el modelo pueden obtenerse 
de los inventarios existentes en las Administraciones de Carreteras. Es asf, como el 
metodo propuesto es de facil aplicacion a parti r de pocos datos asequibles de obtener en 
un puente concreto. 

7. El modelo propuesto es una buena herramienta para ser utilizado a gran escala. Con los 
resultados del modelo, se puede obtener una primera aproximacion de la vulnerabilidad 
de los puentes de una region y, de esta manera los resultados pueden utilizarse para 
identificar los puentes mas necesitados de atencion inmediata. A1 mismo tiempo, se 
pueden estimar que puentes requieren un estudio mas detallado para identificar 
exactamente el tipo de tratamiento que requiere. 

8. Los resultados que se puedan obtener al aplicar el metodo a un area urbana o una region 
determinada, son una buena herramienta para incorporar en los programas de 
planificacion de sistemas de transporte o bien para programas de atencion de emergencia 
en casos de terremoto; ya que el modelo ofrece un estudio preliminar de la vulnerabilidad 
sfsmica de los puentes, definiendo que puentes coirespondes a los mas vulnerables. 
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